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Abstract 


Doppler  shift  measurements  of  the  Na  Dj  absorption  line  have  revealed  solar  oscilla¬ 
tions  in  a  new  regime  of  frequency  and  wavenumber.  Oscillations  of  vertical  velocities 
in  the  temperature  minimum  and  low  chromosphere  of  the  Sun  are  observed  with  fre¬ 
quencies  ranging  up  to  9.5  mHz.  There  is  no  evidence  for  chromospheric  modes  of 
3  minute  period.  This  indicates  that  the  chromosphere  does  not  form  a  good  cavity 
for  acoustic  waves. 

The  fundamental-modes  appear  with  wavenumbers  up  to  5.57  Mm"*  (equivalent 
spherical  harmonic  degree,  3S77).  The  fretiuencies  lie  below  the  predicted  values  at 
wavenumbers  above  I  Mm"*.  The  values  are  in  agreement  with  previous  measure¬ 
ments  that  exist  for  wavenumbers  up  to  2.67  Mm"*. 

Spatial  maps  of  velocity  power  show  that  high  wavenumber  oscillations  are  sup¬ 
pressed  in  active  regions.  The  shape  of  the  power  depression  indicates  that  wave 
motion  is  affected  in  the  layer  of  atmosphere  where  the  nieasurement  is  made.  The 
/-modes  are  suppressed  in  the  same  way  as  p-modes,  indicating  that  the  mechanism 
for  wave  suppression  affects  velocity  fluctuations.  Mode  frequencies  are  not  affected 
by  the  magnetic  fields  by  more  than  50  pHz,  the  precision  of  the  measurement. 
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Chapter  1 
Introduction 


Solar  oscillations  were  first  observed  1)\‘  Leighton  ( 1960).  and  were  explained  by  Ulrich 
(1970)  and  Leibacher  k.  Stein  ( 1971)  as  being  the  superposition  of  many  global  stand¬ 
ing  waves.  They  have  siiljsequently  become  an  important  tool  for  the  determination 
solar  structure.  Measurements  from  helioseismology  include  the  depth  of  the  con¬ 
vection  zone  (Gough  1976).  helium  abundance  in  the  convection  zone  (Christensen- 
Dalsg£uird  &:  Gough  1980),  and  internal  rotation  rate  (Duvall  et  al.  1984;  and  others). 

In  the  past,  observational  limitations  have  restricted  helioseismology  to  long  wave¬ 
length  oscillations  (~  4  Mm).  .New  instruments  and  techniques  are  beginning  to  make 
shorter  wavelengths  accessible  to  study.  This  opens  the  door  to  the  exploration  of 
structure  at  much  smaller  spatial  scales.  Stratification  of  the  atmosphere,  and  vertical 
structure  of  magnetic  regions  are  becoming  active  areas  of  investigation. 

The  term  solar  oscillations  conventionally  refers  to  two  types  of  waves.  Because 
they  exist  in  a  resonant  cavity,  these  waves  are  also  called  modes  (eigemnodes  of 
the  cavity).  Compressional  or  acoustic  waves  (p-modes)  are  longitudinal  waves  with 
pressure  as  the  restoring  force.  They  have  at  least  one  nodal  surface  in  the  radial 
direction.  Gravity  waves  have  gravity  as  the  restoring  force.  A  special  case  of  these 
are  the  surface  waves.  B<‘cause  they  have  no  radial  nodes,  surface  waves  are  called 
fundamental  modes  (/•modes).  These  transverse  waves  propagate  at  the  interface 
between  two  media.  The  remainder  of  the  gravity  waves  are  called  g-taoda.  To  date, 
p-modes  have  not  been  unambiguously  detected,  and  they  will  not  be  considered  hare. 
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1.1  Wave  Properties 

1.1.1  Acoustic  Waves 

Sound  waves  can  propagate  vertically  at  frequencies  greater  than  the  £icoustic  cutoff 
or  critical  cutoff  frequency  u,'c-  At  lower  fretiuency,  a  disturbance  moves  the  whole 
atmosphere  up  and  down  in  phase.  Hence,  waves  do  not  propagate.  At  low  horizontal 
wavenumber  k  ^k  <  w)-  the  acoustic  cutoff  can  be  expressed  in  terms  of  the  adia¬ 
batic  sound  speed  c  and  ( he  density  scale  height  H.  For  an  isothermal  atmosphere, 
the  scale  height  is  given  liy 

H  =  -  (1.1) 

19 

and,  the  acoustic  cutoff,  liy 

^  (1-2) 

where  g  is  the  acceleration  due  to  gravity;  and,  7  is  the  adiabatic  index  (Lamb 
1932).  Taking  c  =  7.1  kut/s.  3  =  278  m/s*,  and  7  =  |  in  the  photosphere,  we  find 
H  =  HO  km  and  uid'lTr  =  5.2  mHz  for  k  4.5  Mm~*. 

In  general,  the  acoustic  cutoff  depends  on  wavenumber.  The  relationship  is  given 
by  (Kumar  &  Goldreich  19S9) 

Figure  1.1  shows  as  a  function  of  k. 

Oscillations  with  wavelength  greater  than  a  few  Mm  propagate  around  the  sun 
without  excessive  dissipation.  .A  periodic  boundary  condition  is,  therefcnre,  imposed 
on  these  waves.  The  wave  equation  is  separable  into  radial  wad  angular  parts  and  has 
eigensolutions  of  the  form. 

♦,.,„(r.0,c>.O  =  Rni{r)Yr{9,<l>)e-*^*  (1.4) 

n  =  l,2,3 -  /  =  0,1,2,...  m  =  1,...,/ 

Short  wavelength  oscillations  are  dissipated  before  traveling  around  the  sun.  More¬ 
over,  at  their  small  spatial  scales,  the  medium  can  be  s^pnudmated  as  plaoe-paralkL 
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Plane  waves,  therefore.  pio\  icle  a  more  natural  basis  set  for  the  horizontal  eigenfunc¬ 
tions. 


(1.5) 


The  horizontal  wavenumber  h  = 
/  as  follows. 


|A'|  can  be  related  to  the  spherical  harmonic  degree 


(1.6) 


Ray  paths  for  a  p-mode  are  diagrammed  in  figure  1.2.  As  waves  travel  downward, 
they  encounter  higher  temperatures  and  an  increasing  sound  speed.  This  causes  them 
to  refract  back  upward  where  they  encounter  increasing  values  of  u?e.  When  the  wave 
can  no  longer  propagate  (^'  <  it  is  reflected  back  down.  The  resonance  condition 
is  met  when  crests  of  downgoing  waves  are  superposed  on  crests  of  refracted  waves 
returning  to  the  surface.  For  waves  trapped  in  an  adiabatic  and  polytropic  layer 
(polytropic  index  s  4),  this  condition  gives  a  dispersion  rdatkm  of  approximate 
form. 
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Figure  1.2:  Propagation  of  an  acoustic  wave  in  a  plane-parallel  atmosphere.  The  wave 
is  refracted  as  it  propagates  downward,  and  reflected  at  the  surface.  The  maximum 
depth  reached  by  the  wav  e  is  called  the  lower  turning  point.  From  Gough  (1983). 


+  n  =  1,2,3,...  (1.7) 

where  ^  is  f  he  acceleration  due  to  gravity  at  the  surface. 


This  dispersion  relation  can  be  measured  directly  from  the  ridges  of  power  in  a 
spectral  decomposition  (a  plot  of  power  vs.  u;  and  k)  of  solar  velocities  or  intensities. 
The  e.xact  relation  depends  on  the  variation  of  sound  speed  with  depth. 

The  lower  turning  point  <1  of  the  waves  is  approximately  given  by. 


3(n  +  2) 


(1.8) 


The  excitation  source  of  the  p-modes  is  thought  to  be  turbulent  convection  at 
the  top  of  the  convection  zone.  The  efficiency  of  acoustic  emission  from  convection 
is  proportional  to  appro.ximately  the  eighth  power  of  the  Mach  number  (Goldreidi 
k  Kumar  1988).  Hence,  only  a  few  of  the  granules  on  the  high  vriocity  tail  of  the 
distribution  are  expected  to  contribute  to  the  acoustic  noise  driving  the  p«modes 
(Brown  1991). 


1.1.2  Surfoce  Waves 

Surface  waves  can  propagate  on  the  interface  between  convectivdy  stable  media  (e.g., 
air  over  watar).  On  the  Sun.  the  transition  region  between  chnmiosidtere  and  corona 
I»ovides  such  a  boundary  (Jones  1969).  The  rkmM  of  Avrett  (1992)  dcKiibee  the 


Figure  1.3:  Temperature  and  density  of  the  solar  atmosphere.  From  Avrett  (1992). 


temperature  and  density  discontinuities  shown  in  figure  1.3.  For  purely  vertical  den¬ 
sity  stratification  of  a  compressible  fluid,  surface  waves  obey  the  dispersion  relation 
(Whitney  1958), 

=  (1.9) 

where  g  is  the  local  acceleration  due  to  gravity. 

The  simple  form  of  this  relation  has  prompted  the  renuurk  that,  “There  is  no 
interesting  physics  in  the  /•modes.’'  It  wasn’t  until  /-modes  were  found  not  to  obey  the 
above  rdation  at  high  wavenumber  that  this  conclusion  was  reexamined.  Libbrecht, 
Woodard  it  Kaufman  (1990)  found  /-mode  frequencies  to  lie  bdow  the  predicted 
values  for  wavenumbers  above  1.4  Mm~*.  At  least  three  possible  explanations  exist, 
though  numerical  predictions  have  not  yet  bera  offered. 

Campbell  &  Roberts  (1989)  and  Evans  Se  Roberts  (1990,  1991)  predict  that  a 
horisoBtal  magnetic  field  in  the  chromosphere  wouhl  increase  the  fireqaeDdes  of  the 
/•modes  by  increasing  the  restoring  force  on  the  surface  waves.  The  s%n  of  the 
frequency  shift  may  reverse  at  high  wavenuiid>er,  however,  in  a  more  realistk  model 
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of  the  chromosphere  with  an  appropriate  magnetic  field  structure  (Roberts  1991). 

Rosenthal  (1992)  finds  that  evanescent  waves  dissipated  in  the  corona  would  cause 
power  to  leak  out  of  the  /-modes.  This  would  lower  the  frequency  of  the  modes  at  a 
given  wavenumber  by  an  amount  dependent  on  the  temperature  difference  between 
chromosphere  and  corona. 

Finally,  structures  in  the  atmosphere,  such  as  flux  tubes  and  convection  cells, 
violate  the  assumption  of  liorizontal  homogeneity.  It  is  not  known  whether  this  should 
increase  or  decrease  the  mode  fre<|uencies  (Gough  1992). 


1.2  High  Frequency  Modes 

For  the  purpose  of  discussion  here,  the  high  frequency  regime  includes  modes  with 
uf  >  iJci  where  tuHz.  This  value  of  ijjf  is  the  photospheric  value  for  low 

wavenumbers.  Because  waves  with  ^  >  u>c  can  propagate,  they  are  not  trapped  by 
the  photosphere.  So,  in  the  model  illustrated  by  figure  1.2,  there  should  be  no  upper 
reflecting  surface  at  frequencies  >  •‘).2  mHz.  Removal  of  the  top  surface  of  the  cavity 
eliminates  the  possibility  of  resonance,  and  hence,  the  mode  structure. 

It  was  a  surprise,  tljercfore.  when  tUe  p-mode  ridges  were  found  to  continue  above 
the  acoustic  cutoff  (Jeffries  et  al.  1988;  Libbrecht  1988;  Duvall  et  al.  1991).  Two  expla¬ 
nations  have  been  put  forward  for  the  origin  of  the  high  frequency  ridges.  Balmforth 
&  Gough  (1990)  propose  that  strong  reflection  from  the  base  of  the  corona  ddines  a 
chromospheric  cavity  which  couples  to  the  subphotospheric  cavity.  Waves  with  fre¬ 
quency  high  enough  to  pass  through  the  photosphere  would  still  be  trapped  by  the 
corona. 

The  presence  of  a  chromospheric  cavity  hu  further  consequences,  however.  The 
chromosphere  is  of  order  one  wavelength  in  thickness  for  5-ininute  modes;  and  waves 
of  appropriate  wavelength  A,  will  resonate  in  this  cavity.  The  p-modes  propagate 
almost  vertically  in  the  atmosphere  (angle  of  incidence  0).  Thus,  A,  =  c/ur, 
independent  of  horizontal  wavenumber.  The  chromospheric  modes  should,  therefore, 
have  a  dispersion  relation. 

a;  =  const.  (f'lO) 
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Figure  1.4:  Mode  frequencies  vs.  wavenumber  with  a  chromospheric  cavity  present. 
From  Ulrich  k  Rhodes  (1977). 

Ulrich  &  Rhodes  (1977)  found  that  the  p>mode  ridges  would  merge  with  a  set  of 
chromospheric  modes  at  about  5.5  mHz.  See  figure  1.4.  The  p-modes  would  continue 
on  above  the  5.5  mHz  band  with  an  "avoided  crossing”  between  the  ridges  above  and 
below  the  chromospheric  resonance. 

Dame,  Gouttebroze  Malherbe  (1984)  clum  to  have  seen  the  chromospheric 
modes  in  Ca  11  K-line  observations.  Th«r  power  spectra  are  very  noisy,  however,  and 
it  is  difficult  to  discern  any  definite  ridges  of  power.  Deming  et  al.  (1986)  find  a  peak 
in  the  one-dimensional  power  spectrum  (power  versus  u;)  at  4.3  mHz.  Its  significance 
is  difficult  to  assess,  though,  without  knowledge  of  the  wavenumber  dependence  of  the 
power.  Woodard  k  Libbrecht  ( 1991 )  do  not  find  this  peak  in  similar  one-dimenstonal 
power  spectra. 

The  other  explanation  for  the  high  frequency  ridges  was  offered  by  Kumar  et  al. 
(1990).  They  propose  that  waves  are  emitted  isotropically  from  a  source  near  the 
of  the  convection  zone.  Waves  traveling  directly  toward  an  overl]ring  observed  surface 
interfere  (in  Fourier  space)  with  waves  refracted  upward  from  below  to  produce  “mock 
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modes.”  It  must  be  stressod  that  this  interfeience  is  not  physical.  In  fact,  the  directly 
emerging  wave  trains  and  refracted  wave  trains  may  arrive  at  the  surface  at  different 
times  and  in  different  places.  The  basis  functions,  used  to  represent  the  waves,  can 
only  distinguish  time,  modulo  'Itt/uj,  and  position,  modulo  2ir/k.  Frequencies  of  the 
mock  modes  can  be  used  to  determine  the  depth  of  the  acoustic  source  layer  (Kumar 
&  Lu  1991). 

The  two  models  for  liigh  frecjuency  waves  can  be  distinguished  by  observations 
that  resolve  the  ridge  structure  in  the  -S.-S  mHz  frequency  range. 


1.3  High  Wavenumber  Modes 

Most  of  the  literature  in  helioseismoiog\-  refers  to  the  (spherical  harmonic)  degree 
/  of  a  mode  rather  than  the  wavenumber  k.  By  convention,  low  degree  refers  to 
0  <  /  <  4,  intermediate  degree,  to  4  <  /  <  100.  and  high  degree  to  100  <  /.  The 
oscillations  presented  in  this  dissertation  have  500  <  /  <  3877.  By  analogy  with 
the  electromagnetic  spect  rum.  they  should,  perhaps,  be  called  very  high  or  ultrahigh 
degree,  but  I  will  stick  to  the  term,  high  wavenumber. 

Previous  observations  of  the  /-modes  exist  up  to  wavenumbers  of  2.67  Mm“' 
(Libbrecht,  Woodard  &  Kaufman  1990).  As  mentioned  above,  the  frequencies  do 
not  agree  with  equation  1.9  at  wavenumbers  above  1.4  Mm~*.  The  variation  of  g 
with  height  accounts  for  only  a  0.5%  modification  to  equation  1.9.  This  is  an  order 
of  magnitude  less  than  the  observed  effect.  It  would  be  desirable  to  confirm  the 
observational  result  and  determine  the  trend  at  higher  wavenumbers. 

The  size  of  structures  that  can  be  resolved  by  oscillations  is  roughly  equal  to 
the  wavelength  of  the  modes.  For  active  regions,  spatial  scales  range  from  200  km 
(flux  tubes)  to  25  Mm  (large  sunspots).  The  corresponding  wavenumbers  range  from 
0.25  to  30  Mm~'.  Although  the  high  end  of  this  range  is  not  attainable  with  any 
current  instrument,  there  is  still  much  room  for  investigation. 

Woods  &  Cram  (1981)  have  reported  a  suppression  of  oscillation  aixq>litude8  in 
plage  as  compared  to  quiet  regions.  Tarbell  et  al.  (1988)  and  Title  tt  a/.  (1992)  find  a 
suppreMion  of  both  convective  motions  and  wave  amplitudes.  The  suppresaicni  of  the 
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oscillations  coincides  with  the  boundary  of  the  photospheric  magnetic  field  to  within 
one  wavelength. 

Braun,  Duvall  LaBonte  (1987,  1988)  and  Braun  et  al.  (1992b)  have  reported 
absorption  of  p-mode  power  liy  sunspots.  The  mechanism  for  this  is  not  understood. 
Resonant  absorption  (Hollweg  1988;  Chitre  Sz  Davila  1991),  scattering  to  short  wave¬ 
length  (Bogdan  Iz  Zweibcl  1987),  and  enhanced  viscous  dissipation  (LaBonte  &  Ryu- 
tova  1992)  are  among  tlie  explanations  offered.  It  is  by  no  means  clear  whether 
only  propagating  acoustic  waves  below  the  photosphere  are  absorbed,  or  whether 
evanescent  waves  are  also  affected.  A  spatial  map  of  oscillation  amplitudes  above 
the  photosphere,  where  the  magnetic  fields  are  different,  would  help  untangle  the  two 
effects. 

Global  oscillations  have  Ijeen  found  to  have  frequency  shifts  correlated  with  solar 
activity  (Palle,  Regulo  Roca  Cortes  1989;  Woodard  et  al.  1991).  This  may  be  due 
to  a  change  in  the  shap*-  of  the  cavity  surface  with  sunspots  present.  Waves  of  a 
given  frequency  may  have  rlifferent  wavenumbers  within  magnetic  and  nonmagnetic 
regions.  Global  measurements  would  then  find  a  weighted  average  of  the  two  disper¬ 
sion  relations.  Compari.son  of  liigh  resolution  observations  in  quiet  sun  and  active 
regions  could  resolve  this  issue  by  direct  measurement  of  dispersion  relations  in  the 
two  types  of  regions. 

1.4  Height  of  Observations 

Oscillations  can  be  observed  as  fluctuations  in  intensity  and  line-of-sight  velocity. 
Both  p-modes  and  /-modes  produce  motions  of  gas  particles,  and  cause  Doppler 
shifts  of  absorption  lines.  Pressure  fluctuations  due  to  p-modes  produce  tempera¬ 
ture/intensity  fluctuations.  The  /-modes  also  produce  intensity  fluctuations,  though 
the  reason  is  not  well  understood.  It  may  be  that  the  <>»  100  km  vertical  displacements 
parcels  of  gas  affect  the  process  of  absorption  line  formation. 

Different  parts  of  an  absorption  line  arc  formed  in  different  layers  of  the  atmo¬ 
sphere.  Hence,  the  height  at  which  the  velocities  and  intensities  are  determined 
depends  on  the  positions  in  the  absorption  line  at  which  measurements  are  made.  In 
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Figure  1.5:  Na  D  line  profiles  on  left  h^ulcl  scale.  Contribution  functions  for  emergent 
intensity  on  right  hand  scale.  Thick  lines  show  the  maximum  of  the  contribution 
function  at  each  wavelength.  Thin  lines  show  1/3  of  the  maximum  of  the  contribution 
function.  From  Schleicher  (1976). 

the  core  of  deep  lines  such  as  Na  0|.  very  little  continuum  light  emerges  from  the 
photosphere.  The  height  of  a  measurement,  therefore,  depends  in  the  contribution 
function  for  the  emergent  line  emission  (Gurtovenko,  Sheminova  L  Sarychev  1991; 
Schleicher  1992).  This  function  is  shown  in  figure  1.5. 

The  original  intention  of  the  observations  described  in  chapter  4  was  to  make  one 
set  of  measurements  in  the  core  and  one  in  the  wings  of  the  Na  Di  line.  Problems 
described  in  section  4.1  resulted  in  a  measurement  weighted  mostly  in  the  core  (600- 
850  km),  but  with  srune  contribution  from  below  the  temperature  miniimim  (150- 
500  km). 


Chapter  2 

Measurement  Techniques 


Producing  power  spectra  from  uncalibrated  images  involves  several  stages  of  data 
analysis  and  reduction.  This  chapter  reviews  some  general  procedures,  st2uidard  and 
nonstandard,  for  measuring  line-of-sight  components  of  velocity  and  magnetic  field, 
and  for  determining  wave  amplitudes  and  phases.  The  details  of  calibration  to  remove 
artifacts  of  the  instrument  and  atmosphere  are  left  to  chapter  5. 

The  line-of-sight  velocity  is  found  from  the  Doppler  shift  of  an  absorption  line. 
The  Doppler  shift,  in  tuin,  is  determined  by  measuring  intensity  as  a  function  of 
wavelength  at  a  series  of  fixed-wavelength  positions  in  the  absorption  line.  The  more 
points  in  the  measurement,  the  better  the  velocity  estimate  will  be.  Unfortunately, 
more  points  require  proportionately  more  time  to  make  the  measurements,  and  geo¬ 
metrically  more  time  to  analyze  them. 

In  this  chapter  are  described  two  new  techniques  for  making  velocity  images  (dopp- 
lergrams)  and  line-center  intensity  mea.surements  using  three-point  and  four-point 
methods.  The  four-point  method  was  developed  for  the  Solar  Oscillations  Investi¬ 
gation  /  Michelson  Doppler  Imager  (SOI/MDI)  experiment  for  the  Solar  and  Helio¬ 
spheric  Observatory  (SOHO)  spacecraft.  It  was  not  applied  to  the  data  used  for  this 
dissertation. 

Images  of  line-of-sight  magnetic  field  strength  (magnetograms)  are  an  important 
tool  in  this  study.  The  tried  and  tested  “two-point  magnetogram”  method  is  described 
below. 
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A  (mA) 

Figure  2.1:  Measurements  of  the  Na  Di  line 

Finally,  the  Fourier  transform  is  essential  to  the  study  of  waves  propagating  in 
a  linear  medium.  Procedures  used  to  extract  information  about  the  oscillations  are 
described  in  the  last  section  of  this  chapter. 

2.1  Three-Point  Measurements  of  Line  Parame¬ 
ters 

The  three-point  method  for  finding  Doppler  shifts  and  line  center  intensities  assumes 
a  parabolic  absorption  line  profile.  This  assumption  is  valid  when  the  three  intoisity 
measurements  are  positioned  relatively  close  to  the  line  core.  Figure  2.1  shows  the 
meMurements  schematically  for  the  Na  Di  (5896  A)  line.  The  Doppler  shift  and  line 
center  intensity  are  the  x  and  y  coordinates  of  the  vertex,  respectivdy. 
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2.1.1  Algebraic  Solution 

Assuming  a  parabolic  line  profile 


\  const. ) 


(2.1) 


we  can  solve  for  the  vertex  (  Aq.  /o)  in  terms  of  three  intensities  /i,  /2,  /a  measured  at 
wavelengths  Ai,A2,A3,  respective!}'. 


Ao  = 

A)  =  Ii  + 


AiCi  +  A3C3 
Cj  +  C3 

(61  Cl  —  ^aCa)^ 


461^3(^3  —  ^i)(C3  —  Cl) 

A3  =  §(A3  +  A2), 

(Si  =  Ai  —  A2.  <^3  =  A3  —  A2, 

<  I  ^  8\{Iz  —  /j).  and  C3  =  —  Ii)- 


(2.2) 

(2.3) 


where  A,  =|(Ai  +  A2), 


2.1.2  Error  Estimates 

The  errors  on  the  measured  quantities  Ix,l2,l3  are  due  to  limitations  of  the  instru¬ 
ment.  The  transmission  of  the  tunable  filter  varies  with  wavelength  as  shown  in 
figure  3.8.  The  phase  of  each  Fourier  component  of  this  modulation  is  independent, 
and  drifts  with  time.  In  addition,  the  wavelengths  tuned  by  the  filter  are  quantized 
in  units  of  the  smallest  step  size  (3.14  mA).  These  sources  of  error  were  used  in  a 
simulation  to  estimate  errors  for  the  velocity  and  intensity. 

The  simulation  assumed  a  parabolic  line  profile  with  parameters  measured  from 
the  data.  A  range  of  velocities  of  ±2000  m/s  with  a  flat  distribution,  and  a  range  of 
line-center  intensities  of  ±8%  with  a  Gaussian  distribution  were  used  in  the  model. 
Intensity  “measurements"  /i,  I2, 13  were  modeled  assuming  a  random  jitter  of  up  to 
±1.57  mA  (flat  distribution)  in  the  location  of  each  measurement,  and  assuming 
phases  for  the  transmission  modulation  which  varied  randomly  between  sets  of  3 
flltergrams,  but  not  within  a  set.  The  velocities  and  intensities  were  then  calculated 
using  equations  2.2  and  2.3.  The  errors  had  approximately  Gaussian  distributions 
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with  standard  deviations. 


(7i,  =  300  m/s,  and  «t/  =  7.7%. 


(2.4) 


2.2  Four-Point  Measurements  of  Line  Parame 
ters 


An  algorithm  was  developed  for  the  SOI/MDI  experiment  to  determine  Doppler  shifts 
using  four  measured  intensities.  The  measurements  are  used  to  solve  a  system  of  si¬ 
multaneous  equations  for  the  four  absorption  line  parameters.  Since  computing  time 
must  be  kept  to  an  absolute  minimum,  the  nonlinear  equations  are  solved  (numeri¬ 
cally)  using  a  look-up  table  approach. 

The  model  included  only  the  Lyot  and  Michelson  filter  elements,  and  the  following 
parameters  were  used: 


•  Lyot  filter:  profile  —  Gaussian 

fwhm  —  500  mA  (see  section  2.2.3) 

•  Michelson  elements  (2):  profile  —  cos*(^'A)  cos*(2fcA) 

fwhm  —  105  mA 

'^here  k  =  ^cos"'  ^2"?  [(1  +  vH)’  +  (^  "  >/i)^]) 

•  Ni  I  (6768  A)  Line:  profile  —  Gaussian 

fwhm  —  130  mA 
depth  —  58% 

•  Measurement  Spacing  —  85  inA  (symmetric  about  nominal  line  center) 

See  figure  2.2  for  illustration. 

Four  measurements  are  made  at  increasingly  longer  wavelengths  across  the  ab¬ 
sorption  line,  to  yield  intensities  Ii  through  I4.  The  ratios,  a  and  are  formed  as 
shown: 


(2.5) 
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Figure  2.2:  Spectral  profiles. 
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/«  _  /- 


(2.6) 


h-h 

The  quantities  a  and  3  are  independent  of  CCD  dark  current  and  gain  variations. 
They  are  also  independent  of  continuum  intensity  and  line-depth  if  a  suitable  choice 
is  made  for  the  Lyot  filter  profile  (see  section  2.2.3). 

We  can  solve  exactly  for  the  Doppler  shift  Aq,  line  depth  d,  and  continuum  inten¬ 
sity  lo  if  we  assume  a  cosine  shaped  absorption  line  with  measurement  spacing  equal 
to  the  hvdtm: 


Solving,  we  find, 


/(A)  =  /o 

(2.7) 

Aq  = 

:  fwhm  •  —  cof* 

(2.8) 

i 

fo  *  2^Ii  +  /s  —  ■*) 

(2.9) 

(2.10) 

where  s  =  yf{Ii  -  /s)*  +  (A  -  /j)’ 


valocify  (m/s) 


Figure  2.3:  Ratio  a  vs.  velocity  for  width  increments:  0,  10,  20,  30%. 


Note  that  a  has  a  one-to-one  correspondence  with  the  ratio  r  = 
which  appears  in  the  SOI/MDI  proposal  (i.e.,  r  =  (q  — 1)/(q-M)).  Hence  determining 
the  velocity  from  a  alone  is  equivalent  to  determining  it  from  tan~‘  r. 

For  a  Gaussian  line-profile, 

/(A)  =  /o  [l  -  (2.11) 

where  «r  =  /wAm/2.36 

The  ratios  a  and  0  are  plotted  as  functions  of  velocity  for  different  line-widths  in 
figures  2.3  and  2.4,  respectively.  It  is  apparent  that  the  velocity  determined  frmn 
a  al<me  would  be  in  error  if  the  calibration  table  were  not  made  for  the  ‘‘right” 
line- width. 

This  caveat  can  be  removed  by  determining  the  vdodty  v  frmn  both  a  and  fix 
V  a:  v(a,  $).  An  algorithm  for  doing  this  is  outlined  below. 


velocity  (m/») 


Figure  2.4:  Ratio  3  vs.  velocity  for  width  increments:  0, 10,  20,  30%. 

2.2.1  Implementing  the  Algorithm 

Look'Up  Table  Approach 

It  is  possible  in  principle  to  compute  velocity  by  using  a  two-dimensional  look-up  table 
in  a  and  0  (i.e.,  v  =  v{q,0)).  This,  howeva,  would  require  an  inordinate  anmunt 
of  memory  ((12  bits  in  a)  x  (9  bits  in  0)  =  2M  words)  to  achKve  an  accuracy  of  a 
few  m/s.  A  better  strategy  is  to  take  advantage  of  the  fact  that  a,  alone,  provides  a 
pretty  good  estimate  of  velocity;  and  write: 

t’(a,  0)  =  Vi(a)  +  v,(q,  0)  (2.12) 

where  va  is  a  correction  to  the  velocity  vt  and,  hence,  does  not  vary  with  a  as  quiddy 
as  does  vt. 

This  approach  can  achieve  an  accuracy  similar  to  the  previous  approach  but  with 
much  less  memory  overhead.  The  table,  vi(a),  requires  a  12  bit  address;  and  the 
table,  V2(a,0),  requires  7  bits  in  a  by  9  bits  in  0  (this  has  not  been  optimiaed);  for 
a  total  of  66k  words. 
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Domain  and  Range  of  a  and  0 

The  ratio  a  has  a  pole  in  the  negative  velocity-domain  (i.e.,  a  red-shifted  line).  To 
avoid  this,  a  test  is  made  to  determine  the  sign  of  the  velocity.  For  positive  velocities: 
a  =  For  negative  velocities:  a  = 

The  velocity  is  negative  if  and  only  if:  Ii  —  J3  >  I4  —  I2  (2.13) 


As  can  be  seen  from  figure  2.3,  a  covers  a  more-or-less  synunetric  r2uige  about 
zero.  Therefore,  it  can  be  scaled-up  to  provide  an  address  for  a  look-up  table  without 
having  to  subtract  an  offset. 

The  function  /?  is  symmetric  in  the  velocity  domain  and  is  continuous  for  sym¬ 
metric  line  profiles.  For  asymmetric  lines  (and  noisy  lines)  0  h2is  a  pole  new  r  =  0. 
In  practice  this  does  not  cause  problems  because,  near  u  =  0,  t;  is  independent  of  0 
(i.e.,  V2{a,0)  »  0  for  small  u).  An  offset  may  have  to  be  subtracted  from  0  when 
scaling  it  for  use  as  an  address. 

Finally,  note  that: 


l2-h  (/l-/3)-(/4-/2) 

h-h  (h-I*) 


(2.14) 


Hence,  both  o  and  0  can  be  obtained  from  the  differences  (/i  —  I3),  (I\  —  I4),  and 

{I4  -  h). 


The  Algorithm 

It  is  assumed  here  that  the  differences:  Ai  =  A  —  A,  A2  =  Ii  —  A,  and  A3  =  I4  — 12 
have  been  accumulated  in  image  buffers. 

For  all  pixels: 


Calculate  0  and  sign  of  velocity: 

numerator  =  Ai  —  A3 

—1.  if  numerator  >  0; 

otherwise. 


stgn 


••I 

_f-l.  if 
“l+l,  o 


0  _  nltm^r4lvr  _ 
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Calculate  a: 


if  sign  >  0: 
otherwise. 


Calculate  velocity: 

V  =  sign  ■  [Vi(C2a)  +  V'2(C3Q,C4 |8)] 


Cl, . . .  ,C4  are  constants;  V''i  and  V2  are  look-up  tables. 

It  should  be  noted  that  noise  in  the  data  may  cause  a  “division  by  zero”  error  in 
the  algorithm. 


Constructing  the  Calibration  Tables 

The  following  procedures  are  used  to  construct  the  calibration  tables  Vi  and  V2. 
Parameters  used  for  this  study  appear  in  square  brackets. 

The  table  Vila)  is  made  for  a  line-prohle  of  nominal  width  [130  mA]  as  follows: 

1.  Calculate  a  for  a  representative  number  [44]  of  discreet  (evenly  spaced) 
velocities  ranging  between  some  number  slightly  less  than  zero  [-40  m/s] 
and  some  number  slightly  larger  than  the  largest  velocity  expected  from 
the  data  [4500  in/s]. 

2.  Invert  the  table  of  o  versus  v  to  get  a  table  of  v  versus  q. 

3.  Interpolate  (cubic  spline)  between  the  endpoints  to  get  V'i(a),  a  table  of 
velocities  for  an  evenly  spaced  abscissa  a. 

The  table  V2(a,l3)  is  made  for  the  same  range  of  velocities  as  is  Vi;  and,  for  a 
given  range  of  widths  [100  to  295  mA].  The  minimum  width  should  not  be  less  the 
measurement  spacing  [S5  mA]. 

There  are  two  “unphysical”  regions  in  the  table  V2.  The  highest  values  of  a  cannot 
occur  with  the  highest  values  of  and  the  lowest  values  of  a  cannot  occur  with  the 
lowest  values  of  These  portions  of  the  table  will  not  normally  be  addressed  by  the 
Doppler  algorithm. 

V3  b  calculated  as  follows: 


CHAPTER  2.  MEASUREMENT  TECHNIQUES 


20 


1.  Calculate  a  and  ^  for  a  representative  grid  of  velocities  [22]  and  widths 
[14]. 

2.  Interpolate  (spline)  to  get  a  much  finer  grid  [512  in  velocity  by  256  in 
width]. 

3.  For  each  “real"  velocity  v  [512  of  them]  calculate  the  velocity  correction 
t’2: 

v,  =  y  -  Vi(o) 
using  the  table  HCo)  computed  above. 

4.  Sort  the  ordered  triplets  (w2^  Q,  S)  into  a  2-dimensional  table  ^(a,  0),  keep¬ 
ing  track  of  how  many  times  each  bin  gets  filled. 

5.  Fi.x  uncooperative  bins: 

a)  If  a  bin  accumulated  more  than  one  value  of  normalize  to  get 
the  “average”  Uj  for  that  (o,  /?). 

b)  If  a  bin  in  a  “physical”  region  did  not  get  filled  [there  were  6  of 
them],  interpolate  [linearly]  using  adjacent  bins. 

c)  Fill  “uiiphysical”  bins  with  the  value  of  vj  from  the  nearest  “physi¬ 
cal”  bin.  This  is  done  because  noise  in  a  measurement  may  produce 
an  ordered  pair  (a. /d)  which  does  not  lie  quite  where  it  should  in 
the  table. 

Tables  Vj  and  Vj  were  made  in  units  of  ro/s,  and  rounded  to  the  nearest  integer. 

2.2.2  Spectral  Simulation 

The  instrument  and  Doppler  algorithm  were  simulated  in  order  to  determine  the 
magnitude  of  statistical  and  systematic  errors.  The  simulation  program  stored  the 
line  and  transmission  profiles  in  single-precision  512-element  tables  covering  the  range 
±800  mA  about  nominal  line-center.  The  intrasities,  /i, . . . , 7^,  were  calculated  for  a 
discrete  set  of  velocities  as  follows: 

1.  Offset  the  index  of  the  line- profile  Ud>le  to  simulate  a  Doppler  shift. 
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2.  Take  the  product  of  the  transmission- profile  and  the  offset  line- profile. 

3.  Integrate  [sum]  over  wavelength  to  get  “measured”  intensity. 

The  Riemann  interval  and  limits  of  integration  were  varied  to  verify  that  they 
did  not  affect  t  he  results. 

The  error  estimates  are  described  in  the  next  two  sections. 


2.2.3  Sources  of  Systematic  Error 

Line-Depth  Variations 

The  ratio  a  will  be  sensitive  to  variations  in  line  depth  if  measurements  /i  and  I3  (or 
I4  and  I2)  do  not  integrate  the  same  amount  of  continuum. 

Figure  2.5  shows  the  transmission  through  the  Lyot  and  Michelson  filters  as  a 
function  of  wavelength  for  the  first  and  third  measurements.  The  ratio  of  integrals  of 
these  curves  is  plotted  as  a  function  of  Lyot  filter  fwhm  (assuming  a  Gaussian  Lyot- 
profile)  in  figure  2.6.  The  two  integrals  are  equal  when  the  Lyot  fwhm  is  625  mA. 
They  are  different  by  2%  when  the  fwhm  is  500  mA. 

This  2%  difference  causes  a  velocity  error  £„which  depends  on  line-depth  as  shown 
in  figure  2.7.  If  the  Lyot  fwhm  were  400  mA,  e„  would  be  larger  by  a  factor  of  4. 

Lyot  Filter  Drift 

If  the  Lyot  filter  profile  drifts  in  wavelength  space,  the  four  measured  intensities  may 
not  integrate  the  same  amount  of  continuum.  This  would  give  erroneous  values  for  a 
and  0  and,  hence,  v.  To  eliminate  this  source  of  error  the  transmission  profile  t/,i,t(A) 
should  satisfy: 

J  =  constant  i  =  l,...,4  (2.15) 

for  some  small  range  over  which  the  Lyot  filter  might  drift;  t\^ici^{X)  is  the  trans¬ 
mission  of  the  Michelson  filter  at  wavdength  setting  i ; 

Fcm:  a  Gaussian  Lyot-profile,  it  turns  out,  this  condition  is  more  <a  less  satisfied 
for  the  same  fwhm  that  minimizes  due  to  line-depth  variations.  The  reascm  can  be 
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-800  -400  0  400  800 

Wavelength  (mA) 

Figure  2.5:  Transmission  product  for  Lyot  and  Michelson  filters  vs.  wavelength.  U 
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Figure  2.7:  Systematic  error  c„  vs.  velocity  for  line-depth  increments 
-20,-10,0,10,20%. 

seen  from  figure  2.5.  For  transmission  profile,  ti,  if  the  Lyot  profile  shifts  to  the  left 
(right),  any  increase  (decre2ise)  in  the  area  under  the  main-lobe  is  compensated  by  a 
decrease  (increase)  in  the  area  under  the  largest  side-lobe.  Transmission  profile  tz  is 
unchanged  to  first  order  since  it  lies  close  to  the  maximum  of  the  Lyot-prolile. 
Figure  2.8  shows  for  S  mA  offsets  of  the  Lyot  filter. 

Line-Width  Variations 

As  stated  above,  if  a  alone  is  used  to  compute  v,  an  error  will  be  introduced  which 
depends  on  the  line- width.  Figure  2.9  shows  the  error  incurred  if  the  line- width  is  10, 
20,  and  30%  larger  than  the  nominal  width  for  which  the  calibraticm  table  was  made. 
If  both  a  and  $  are  used,  the  error  is  reduced  to  a  few  m/s  as  shown  in  figure  2.10. 


error  V(mdi) 


CHAPTER  2.  MEASUREMENT  TECHNIQUES  25 


velocity  (m/s) 

Figure  2.10:  Systematic  error  €„  vs.  velocity  for  line- width  increments:  0, 10,  20,  30% 
(a,  0  method). 

Line- Asymmetries 

The  effect  of  line-asymmetry  was  explored  by  using  a  line-profile  with  a  Gaussian- 
shaped  bisector  as  shown  in  figure  2.11.  It  was  found  that  the  depth  at  which  the 
velocity  is  measured  depends  on  the  velocity,  but  is  independent  of  bisector  shape. 
The  velocity  dependence  is  plotted  for  different  bisector  shapes  in  figure  2.12  (a-only 
method),  and  figure  2.12  (a,  3  method).  Note  that  ty  does  depend  on  the  precise 
shape  of  the  bisector,  but  its  magnitude  cannot  be  determined  until  a  measurement 
is  made  of  the  bisector  shape  at  the  spatial  resolution  of  the  MDI  instrument. 

2.2.4  Statistical  Error 

Statistical  errors  were  found  by  a  standard  propagation  of  error  calculation.  Deriva¬ 
tives  of  the  functions  vt(o)  and  1/2(0^ 3)  were  approximated  by  finite  differences. 
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velocity  (m/s) 

Figure  2.13:  Depth  of  measured  velocity  vs.  velocity  for  bisector  shape  hwhm: 
—7.5, . . . ,  7.5  mA  (a,  i3  method) 


For  a  sinusoidal  line-prohle  with  delta-function  transmission-profiles,  the  statisti¬ 
cal  error  <r„  can  be  written  down  in  a  form  like  that  shown  in  the  SOI/MDI  proposal. 


2cs  1  y/2^ 

d 


(2.16) 


where,  s  is  the  measurement  spacing;  d,  the  line  depth;  and,  Iq,  the  number  of  photons 
collected  at  continuum  intensity  (half  the  well  depth  for  the  CCD). 

Numerically,  ety  from  this  calculation  is  16  m/s  (with  d  =  .58,  /q  =  10’).  The  finite 
spectral  resolution,  howe\'er,  tends  to  wash  out  the  signal.  For  Gaussian  transmission- 
profiles  (Judim  ss  100  mA),  is  25  m/s.  For  the  transmission-profiles  in  figure  2.5, 
<T„  is  30  m/s  as  shown  in  figure  2.14  (a-only  method)  and  figure  2.15  (a,  0  method). 
The  scatter  in  figure  2.15  is  due  to  the  method  by  which  the  calibration  table  ^3(0, 0) 
is  calculated. 

The  sudden  increase  of  <7^  for  velocities  above  3500  m/s  is  due  to  the  increasing 
dependence  of  i;  on  d  together  with  the  fact  that  0  is  only  allotted  9  bits  in  the 
calibration-table  address.  Allotting  more  bits  to  0  would  proportbnately  reduce  this 
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behavior. 


2.2.5  Discussion  of  the  Method 

The  studies  indicate  that  velocity  errors  due  to  line-depth  variations  and  Lyot  filter 
drift  can  be  minimized  by  proper  design  of  the  Lyot  filter  profile.  In  principle  it 
should  be  possible  to  eliminate  the  errors  entirely,  though  only  a  crude  optimization 
has  been  done  here. 

Line- width  variations  can  produce  velocity  errors  of  up  to  a  few  hundred  m/s 
when  the  velocity  is  calculated  from  a  single  rational  function  of  the  data  (a-only 
method).  The  error  is  reduced  to  a  few  m/s  by  using  two  rational  functions  (o, 
method).  The  more  accurate  computation  requires  an  extra  image  buffer  and  entails 
doing  2.2  times  as  many  operations. 

A  line-asymmetry  will  affect  the  measured  velocity  because  the  line- bisector  is 
sampled  at  different  depths  across  the  solar  disk.  This  happens  because  the  average 
velocity  is  a  function  of  disk  position,  and  the  data  tend  to  sample  more  of  the 
line-wings  for  large  Doppler  shifts. 

Data  taken  at  La  Palma  in  19SS  show  that  the  averaging  of  subresolution  velocity 
elements  produces  a  “noise”  level  of  ±20  m/s.  Hence,  it  seems  reasonable  to  set 
this  as  an  upper  limit  for  systematic  errors. 

The  statistical  error  for  a  four  point  measurement  is  ~  ±30  m/s  with  a  CCD 
well-depth  of  2  x  10*  photons.  The  error  becomes  significantly  larger  at  v  ~  3500 
m/s  for  the  a,  0  method.  This  can  be  rectified  by  use  of  a  larger  calibration  table. 
The  errors  for  2  pass  (8  point)  and  4  pass  (16  point)  measurements  are  21  m/s  and 
15  m/s,  respectively. 

In  order  to  minimize  the  statistical  error,  sidelobes  should  be  kept  to  a  minimum. 
Hence  the  optimum  Lyot  filter-profile  would  have  a  flat  peak  with  wings  just  large 
enough  to  ensure  integration  of  equal  amounts  of  continuum  at  each  wavelength 
setting. 

This  recommendation  for  the  Lyot  filter  width  has  been  used  in  the  design  of  the 
MDI  instrument.  The  choice  of  velocity  calculation  method  is  still  pending  furth«' 
study.  It  is  not  yet  known  whether  line-width  or  line-asynunetry  variations  will 
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contribute  more  to  measurement  errors.  Only  if  line-asymmetry  variations  are  small 
can  the  a,  0  method  used. 

2.3  Magnetograms 

2.3.1  Inverse  Zeeman  Effect 

Spectral  lines  of  atoms  in  a  magnetic  field  are  split  by  the  2^man  effect.  For  a  field 
component  pointing  toward  the  observer,  an  emission  line  of  wavelength  A©  is  split 
into  two  components  (<t  components)  with  opposite  circular  polarizations.  The  left 
and  right  circular  polarizations  (LCP  and  RCP)  are  displaced  from  A©  by  — Aa  and 
-I-Aa,  respectively. 

Aa  =  -^A©»(5'A/'  -  g"M'')  (2.17) 

iirnieC* 

where  g'  and  g"  are  the  Lande  factors,  and  Af'  and  M"  are  the  magnetic  quantum 
numbers  for  the  initial  and  final  states  (Bray  k  Loughhead  1964).  The  magnetic 
quantum  numbers  are  subject  to  the  selection  rule  M"  —  M'  =  ±1.  For  the  Fe  I 
(6302.508  A)  line, 

^  =  47.4  mA/kgauss.  (2.18) 

"II 

For  absorption  lines,  the  LCP  and  RCP  components  will  be  preferentiaUy  absorbed 
at  A©  —  Aa  and  A©  -l-  Aa,  respectively.  This  is  termed  the  inverse  Zeeman  effect. 

2.3.2  Longitudinal  Field  Measurement 

From  figure  2.16  we  can  write  down  the  measured  LCP  and  RCP  intensities. 

/t  =  /o+^AA  (2.19) 

/r  =  /o-Kaa  (2.20) 

where  is  assumed  to  be  constant  over  the  range  A.,  ±  Aa- 
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We  can  now  define  the  magnetogram  signal  M, 


M  = 


h  ~  Ift 
h  +  Ir 


/o 


Combining  this  with  equation  2.18,  we  have. 


(2.21) 


H|l  =  cM  (2.22) 

where  the  constant  of  proportionality  c  depends  on  the  profile  of  the  absorption 
line  and  of  the  filter  through  which  the  intensities  are  measured.  The  constant  is 
independent  of  continuum  intensity.  For  the  filter  pariuneters  in  table  4.1, 

c  «  10^  gauss.  (2.23) 


2.4  IVansform  Methods 

The  goal  of  the  data  analysis  is  to  measure  the  dispersion  rdation  and  amplitude 
of  waves.  The  natural  method  for  plane  wave  decompositkm  in  a  linear  medium  is 
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the  Fourier  transform.  Consider  a  monochromatic  plane  wave  p  with  wave  vector 
ko  =  {kox,koy)  and  frequency  u;o: 

p{x,  t)  =  cos{kQ  ■  X  -  u!ot)  (2-24) 

The  Fourier  transform  of  p  is  defined  as, 

/oo  too  too  ^  - 

I  I  cos(^  •  X  —  dt 

'00  «/-oo  •/— oo 

Evaluating  the  integral, 

=  \^Hk  -  ko)  6(u:  -  wo)  +  ^S^{k  +  1^)  6(w  +  wo) 

The  two  delta  functions,  l\  ing  opposite  one-another  in  (^,  w)  space,  illustrate  a  general 
property  of  Fourier  transforms.  For  real  valued  functions,  the  (complex)  Fourier 
transform  will  have  a  factor  of  2  redundancy  in  information.  In  fact,  we  only  need 
to  look  at  a  half-space  (containing  the  origin)  to  obtain  all  the  information  about 
the  original  function.  It  is  conceptually  easiest  to  eliminate  the  negative  frequency 
half-space,  and  consider  only: 


(2.25) 

(2.26) 


— oo  <  kj.  <  00.  — oo  <  ky  <  oo,  0<w<oo 

It  is  important  to  account  for  relative  motion  of  the  observers  field  of  view  with 
respect  to  the  observed  medium.  If  the  medium  is  moving  at  constant  velocity  v  with 
respect  to  the  observer,  we  have  x  —*  x  —  vt.  This  produces  a  measured  frequency 
shift  u;  — m;  4-  for  the  plane  wave,  where 

X,  =  kv.  (2.27) 

A  final  and  extensively  used  property  of  Fourier  transforms  is  the  convolution 
theorem.  It  states  that  the  Fourier  transform  of  the  convolution  of  two  functions 
/(x)  and  ^(x)  is  equd  to  the  product  of  the  Fourier  transforms  of  the  functions. 

-()<!(}  XWi))  (2.28) 
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2.4.1  Practical  Considerations 

The  discrete  Fourier  translbrni.  or  Fourier  series,  can  be  viewed  as  an  extreme  example 
of  a  Fourier  transform  (Bracewell  1986).  The  series  representation  f*  is  useful  for 
discrete,  evenly  spaced  data  /„  with  periodic  boundaries.  It  is  defined  by, 

F,  =  1  X:  (2.29) 

■'''  n=l 

The  inverse  transform  is  defined  by, 

/„  =  ,F-‘{F,.}  =  1  XI  (2.30) 

^  k=i 

For  measurements  made  at  an  interval  6,  the  maximum  Fourier  frequency  that  cam 
be  represented  is  the  Ny(|uisl  freciuency. 

/nx,  =  ^  (2.31) 

The  discrete  transform  has  the  advantage  that  it  can  be  evaluated  using  a  very 
efficient  algorithm,  the  Fast  Fourier  Transform  (FFT)  algorithm.  The  data  from  the 
series  of  Doppler  and  intensity  images  is  indeed  binned  discretely  on  am  evenly  spaced 
grid  in  the  spatial  dimensions  (by  pixel)  and  temporad  dimension  (by  imatge  number). 
The  condition  of  periodic  boundaries  must  be  enforced,  however,  by  apodizing  (see 
below). 

There  are  several  algorithms  available  for  discrete  Fourier  amd  rdated  tramsforms. 
For  handling  real  valued  data,  it  is  more  efficient  (by  a  factor  of  2)  in  computing 
time  and  memory  usage  to  construct  the  Fourier  tramsform  from  a  Fast  Sine-Cosine 
Transform.  For  doing  convolutions  to  smooth  amd  filter  data,  it  is  most  efficient  amd 
conceptuadly  simplest  to  use  the  Faist  Hartley  Tramsform.  The  abbreviation  FFT  is 
used  generically  hereafter  to  refer  to  the  Sine-Cosine  amd  Ha^tl^  tramsfawms. 

The  FFT  adgorithms  used  were  designed  for  data  series  oontauning  N  points,  where 
N  is  even  and  factorable  into  low  primes.  The  spatiad  dimensions  were  512  pixeb  in 
both  X  and  y  (N,  =:  N^^z  >*).  The  temporad  dimension  was  480  images  (^i  =  2‘'3-5). 
This  number  was  chosen  liecause  it  entiuied  the  largest  ‘‘nice”  number  of  filtergrams 
that  would  fit  on  a  single  data  tape. 
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2.4.2  Apodization 

To  approximate  periodic  boundaries,  the  data  is  multiplied  in  each  dimension  by  a 
function  which  tapers  to  zero  at  the  edges.  The  exact  form  of  the  function  must 
be  chosen  very  carefully  to  avoid  disturbing  the  Fourier  spectrum.  The  convolution 
theorem  tells  us  that  multiplying  the  data  by  an  apodizing  function  is  equivalent  to 
convolving  the  Fourier  spectrum  with  the  Fourier  transform  of  the  apodizing  function. 
Therefore,  the  apodizing  function  should  be  constructed  such  that  its  transform  has 
a  narrow  bandwidth  and  very  small  sidelobes. 

Several  standard  apodizing  functions  are  listed  in  Harris  (1978).  The  3-tenn 
Blackman- Harris  function  was  chosen.  Its  discrete  Fourier  transform  has  a  3  dB 
width  of  1.56  bins,  and  largest  sidelobes  of  —61  dB.  The  windowing  function  for  N 
bins  is  given  by 

H  n  =  «o  +  +  “2  (17”)  (2.32) 

for  n  =  -fl, . . . ,  ^  -  1 

where  oo  =  0.44959 
«i  =0.49364 
fl2  =0.05677 

When  filtering  data  as  done  in  chapter  7,  it  is  not,  in  principle,  necessary  to 
apodize  the  data.  For  the  La  Palma  data  sets,  however,  the  image  rotates  during 
the  observations.  Removal  of  this  rotation  from  the  images  leaves  the  corners  of  the 
images  without  real  information.  To  minimize  round  off  errors  from  the  comers,  the 
data  were  masked  before  filtering.  The  mask  used  had  a  cosine  taper  (to  zero)  at  the 
edges  of  a  flat  central  disk: 


1  ifr„n<iV-32 

0  if  r,„n  >  N 

cos  (^(r„,„  -  N  +  32))  otherwise 


(2.33) 


where  r,n„ 


=  y/rrir+T?, 

=  . f-1 


m,n 
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2.4.3  Normalization 

The  FFT's  are  normalized  such  that  the  spectral  power  of  a  sine  wave  is  equal  to  the 
square  of  the  wave-height  before  apodization. 

f-i 

E  H'Vn  ■  5  sin  kon}\^  =  1  (2.34) 

The  exact  normalization  is  a  weak  function  of  k  and  is  different  for  the  real  and 
imaginary  parts  of  the  transform.  However,  it  can  be  approximated  by  a  constant. 

|jF{/„}p  =  .41:F{VF„/„}p  (2.35) 

(2.36) 

This  approximation  is  valid  for  wavelengths  that  are  short  compared  to  the  scale 
of  variation  of  the  apodizing  function  as  is  the  case  for  all  wavelengths  and 
frequencies  of  interest. 

There  is  an  additional  wavenumber-dependent  normalization  factor  that  corrects 
for  the  attenuation  of  spatial  frecjuencies  by  the  atmosphere  and  optics.  This  correc¬ 
tion  is  described  in  sections  3.4  and  5.4.1. 

2.4.4  Statistical  Error  Estimates 

An  error  estimate  can  be  obtained  for  the  spectral  power  |F*|*  by  a  standard  propa¬ 
gation  of  error  calculation.  If  we  assume  a  constant  error  oy  on  the  data  /„,  we  can 
use  equation  2.29  to  find  the  error  on  the  power 

■'If.P  =  (2.37) 

In  three  dimensions,  this  becomes; 


where  .4  = 


N 

T 


(2.38) 


Chapter  3 

Telescope  and  Instrumentation 


3.1  Telescope 

The  50  cm  Vacuum  Solar  Telescope  at  the  Swedish  Solar  Observatory  is  shown 
schematically  in  figure  3.1.  photograph  of  the  tower  at  Roque  De  Los  Mucha- 
chos  on  the  island  of  Sail  Miguel  De  La  Palma  (Canary  Islands,  Spain)  appears  in 
figure  3.2.  The  telescope  is  a  refractor  with  an  /  number  of  45. 

The  image  forming  lens  is  a  50  cm  achromatic  doublet.  Two  flat  mirrors  Ml  and 
M2  in  the  altazimuth  turret  guide  the  beam  down  an  evacuated  tube.  The  first  lens 
of  the  achromat  forms  the  vacuum  window.  At  the  bottom,  the  flat  mirror  M3  directs 
the  beam  through  a  vacuum  window  to  the  (Ursies)  observing  room.  Extreme  care 
was  taken  by  the  designers  of  the  optics  to  ensure  the  telescope  would  perform  at  the 
diffraction  limit  (Scharmer  tt  al.  1985). 


3.2  Postfocus  Instrumentation 

3.2.1  Tunable  Filter 

Lyot  Filter  Operation 

A  Lyot  filter  is  constructed  from  a  series  of  simple  birefringent  elements.  A  represen¬ 
tation  of  one  such  element  is  shown  in  figure  3.3.  An  incident  pulse  of  light  is  polarised 
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Figure  3.1:  Schematic  diagram  of  the  Swedish  50  cm  Vacuum  Solar  Telescope.  From 
Scharmer  et  ai  (1985). 

at  45”  to  the  optical  axes  of  a  birefringent  crystal.  Fast  and  slow  pulses  emerge  from 
the  crystal,  and  pass  through  an  exit  polarizer.  The  output  of  the  element  is,  .hus, 
a  pair  of  identically  polarized,  equal  amplitude  pulses,  with  a  temporal  separation 
proportional  to  the  crystal  thickness.  The  spectral  response  of  such  an  element  is 
given  by  the  Fourier  transform  of  its  impulse  response.  The  Fourier  transform  of  a 
pair  of  delta  functions  is  a  cosine. 

A  Lyot  filter  consists  of  a  cascade  of  elements,  each  double  the  length  of  the 
previous  one  (see  figure  3.4).  The  transmission  of  each  element  as  a  function  of 
frequency  is  shown  in  figure  3.5.  The  longest  element  has  the  most  closdy  spaced 
maxima.  The  next  longest  has  minima  that  fall  on  every  second  maximum  of  the 
longest  element,  and  so  forth.  The  net  effect  is  to  produce  a  transmission  profile 
with  a  mainlobe  width  approximately  equal  to  that  for  the  longest  element;  and  a 
free  spectral  range,  or  distance  between  peaks,  equal  to  the  period  for  the  sh(»test 
dement. 


Figure  3.2:  Swedish  Solar  Observatory,  La  Palma,  Spain. 
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Figure  3.3:  Pulse  propagation  through  a  simple  element.  From  Title  k  Rosenberg 


Figure  3.5:  Transmission  vs.  frequency  for  the  components  of  a  Lyot  filter  (A,  B, 
C,  D)  and  the  transmission  through  successive  stages  of  the  cascade  (E,  F,  G).  FVom 
Title  Sc  Rosenberg  ( 1981 ). 


CHAPTER  3.  TELESCOPE  AND  INSTRUMENTATION 


41 


Figure  3.6:  Birefringent  tuning  element.  From  Title  k  Rosenberg  (1981). 
Tuning  a  Birefringent  Element 

The  bandpass  of  a  birefringent  filter  can  be  tuned  over  the  full  free  spectral  range 
by  individually  tuning  the  component  elements.  The  best  practical  method  for  doing 
this  is  illustrated  in  figure  3.6  (Title  i:  Rosenberg  1981).  As  in  figure  3.3,  a  fast 
and  slow  pulse  emerge  from  the  crystal  parallel  to  the  fast  and  slow  axes.  The  delay 
between  the  pulses  is  inversely  proportional  to  wavelength.  Reconstituting  the  two 
polarizations  gives  elliptically  polarized  light.  The  axes  of  the  ellipse  are  horizontal 
and  vertical,  and  the  ratio  of  the  axes,  the  ellipticity,  is  a  function  of  wavelength. 

A  quarter  wave  plate  transforms  the  polarization  from  elliptical  to  linear,  and  the 
orientation  of  the  resulting  polarization  vector  is  dependent  on  ellipticity.  Finally, 
a  rotatable  half  wave  plate  and  fixed  exit  polarizer  are  used  to  select  the  linear 
polarization,  and  hence,  the  wavelength. 
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The  Solar  Optical  Universal  Polarimeter 

The  Solar  Optical  Universal  Polarimeter  (SOUP),  developed  by  LPARL,  employs  a 
nine  element  tunable  birelVingent  filter.  It  is  similar  in  design  to  a  Lyot  filter,  but  each 
birefringent  element  contains  a  partial  polarizer  to  apodize  its  transfer  function.  The 
result  is  a  filter  with  smaller  sidelobes  than  a  Lyot  filter.  The  SOUP  instrument  was 
used  on  the  Spacelab  2  shuttle  flight,  and  the  filter  has  become  the  evaluation  model 
for  the  Orbiting  Solar  Laboratory  —  Coordinated  Instrument  Package  (OSL-CIP). 
Software-selectable  blocking  filters,  each  6-S  A  wide,  are  centered  on  solar  absorption 
lines.  Waveplates  can  also  be  selected  to  measure  left  and  right  circular  polarizations 
as  well  as  linear  polarizations  spaced  every  45°. 

At  a  wavelength  of  .\.  the  free  spectral  range  of  the  tunable  filter  is  14  A, 
and  the  fwhm  of  the  passbancl  is  62  m.4.  A  "wide  mode”  is  also  available  by  rotating  a 
half  wave  plate  within  the  longest  element.  This  effectively  shortens  the  element  by  a 
factor  of  3,  and  increases  the  bandwidth  to  103  mA.  The  stepping  motors  that  rotate 
the  tuning  half  wave  plates  i^evmit  tuning  in  3.14  mA  steps.  Tuning  is  repeatable  to 
within  half  a  step. 

Temperature  gradients  and  imperfections  in  the  filter  components  produce  small, 
but  measurable,  variations  in  filter  characteristics  over  the  field  of  view.  The  center 
of  the  passband  varies  b,\-  up  to  10  niA  and  the  fwhm,  by  up  to  5  mA  over  the  field  of 
view.  In  addition,  the  total  amount  of  light  transmitted  (integrated  over  wavelength) 
is  modulated  as  the  filter  is  tuned.  (See  figure  3.8.)  This  modulation  accounts  for 
most  of  the  error  in  the  \elocity  and  intensity  measurements. 

3.2.2  CCD  Camera 

The  CCD  (charge  coupled  device)  camera  uses  a  Texas  Instruments  1024  x  1024  pixel, 
virtual  phase  detector.  It  is  also  an  evaluation  model  for  OSL-CIP.  Subareas  of  the 
full  detector  can  be  read  out  to  reduce  the  time  required  to  write  data  to  tape.  In 
addition,  2  x  2  groups  of  i)i.\els  can  be  summed  to  decrease  exposure  time  by  a  factor 
of  4.  This  2  X  2  summing  morle  was  used  for  the  majority  of  observations  listed  in 
chapter  4. 
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Figure  3.7:  Transmission  profile  of  the  SOUP  tunable  filter  vs.  wavelength.  Trans¬ 
mission  was  measured  using  a  spectrometer,  with  the  filter  tuned  to  6328  A. 

The  detector  geometry  consists  of  a  rectangular  grid  of  pixels.  Each  pixel  is 
efficient  at  collecting  light  up  to  the  boundary  with  its  adjacent  pixeb.  Thus,  the  full 
area  of  the  CCD  array  is  sensitive  to  light. 

The  charge  accumulated  on  a  pixel  is  proportional  to  the  integrated  light  falling 
on  that  pixel.  The  signal  also  integrates  a  current  present  in  the  absence  of  light, 
the  ‘*dark  current.”  The  dark  current  depends  on  device  temperature.  To  keep  this 
effect  constant,  and  to  a  minimum,  the  CCD  is  kept  in  a  vessel  at  low  pressure,  and 
is  cooled  thermoelectrically  to  —20®  C. 

3.2.3  Image  Stabilization  System 

Compensation  is  made  for  small  image  displacements,  caused  by  seeing  and  inaccurate 
telescope  tracking,  by  an  active  tilt  mirror  in  the  beam  path.  The  mirror  is  moved 
in  accordance  with  signals  from  a  sunspot  tracker  which  senses  image  motion.  The 
active  mirror  is  mounted  on  three  piezoelectric  itctuators  which  can  tilt  the  mirror  on 
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Figure  3.8:  Relative  transmission  of  the  SOUP  tunable  filter  vs.  tuned  wavelength. 
Transmission  was  measured  using  a  continuum  illumination  source. 

two  axes.  The  system  has  a  bandwidth  of  65  Hz. 

The  spot  tracker  is  an  array  of  four  photodiodes  (a  quad-cell).  A  sunspot  image 
centered  on  the  quad-cell  will  produce  equal  signals  from  each  diode.  If  the  image 
moves  toward  one  of  the  four  quadrants,  however,  that  quadrant  will  experience  less 
illumination  and  produce  a  smaller  signal.  The  differences  of  signals  from  the  four 
quadrants  are  used  to  determine  compensatory  signals  for  the  active  mirror.  The 
integrated  compensatory  signals  are  used  a  feedback  to  the  telescope  tracking 
system. 

3.3  Optical  Layout 

The  optical  bench  was  set  up  in  the  Ursies  room  as  shown  in  the  photograph  in 
figure  3.9.  The  schematic  layout  of  components  is  shown  in  figure  3.10.  The  beam 
from  the  telescope  strikes  the  active  mirror,  and  is  deflected  through  a  field  stop 
slightly  larger  than  the  active  area  of  the  CCD.  The  heat  stop,  used  to  filter  out 


CHAPTER  3.  TELESCOPE  AND  INSTRUMENTATION 


45 


infrared  wavelengths,  was  removed  from  the  light  path  because  it  produced  fringing 
patterns. 

Wavelengths  of  5000-  7000  A.  are  reflected  by  the  dichroic  beam  splitter  to  the 
tunable  filter  path.  Infrared  wavelengths  are  not  reflected.  Other  wavelengths  are 
reimaged  through  the  Ca  II  K-line  filter  and  the  4308  A  bandpass  filter  to  the  K-line 
and  “white  light”  cameras.  respectivel\\  The  white  light  images  were  used  to  monitor 
the  seeing  conditions  and  the  focus  in  real-time.  The  K-line  camera  was  not  used  for 
these  observations. 

After  reflection  from  the  dichroic  beam  splitter,  the  beam  is  collimated  by  a  lens 
on  a  motor  driven  stage.  The  magnification  of  the  system  is  chosen  by  substituting 
lenses  of  different  focal  h'ngth.  The  stage  is  moved  automatically  to  compensate  for 
the  chromatic  dispersion  of  the  telescope,  and  the  whole  assembly,  moved  manually  to 
adjust  the  focus.  The  beam  pa.sses  through  an  aperture  stop  and  is  reimaged  through 
a  set  of  selectable  waveplates  and  blocking  filters,  and  the  tunable  filter.  For  solar 
observations,  the  beam  passes  through  a  sector  shutter  to  the  CCD  camera.  To  check 
the  position  of  a  spectral  line,  the  flip  mirror  was  inserted  into  the  beam  to  deflect  it 
to  the  photometer. 

Part  of  the  beam  is  di\  erted.  Ijefore  the  polarizer  wheel,  to  the  spot  tracker.  There 
is  a  neutral  den.sity  filter  and  adjustable  zoom  lens  system  in  this  path  to  optimize 
the  intensity  and  magnifi<  ation  of  the  image  for  the  spot  tracker. 

The  motor  driven  stage,  polarizer  wheel,  blocking  filter  wheel,  and  tunable  filter 
were  controlled  by  a  (PC'  compatible)  COMPAQ  computer.  The  shutter  and  CCD 
camera  were  controlled  by  a  \’.\X  Station  3200.  The  system  was  free-running,  with  no 
timebase  to  control  the  data  ac(|uisition.  The  timing  of  the  images,  therefore,  varied 
as  shown  in  figure  3.11.  Tlie  time  mark  given  to  each  image  came  from  a  crystal 
controlled  clock,  and  is  a<  curate  to  ~  I  ms. 
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Figure  ;{.9:  Optical  bench  setup. 
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OPTICAL  LAYOUT  AT  SWEDISH  SOLAR  OBSERVATORY 
LA  PALMA,  SUMMER  1990 


Figure  3.10: 


tima  (min) 


Figure  3.11:  Timing  offset  for  images  vs.  time  from  start  of  run.  The  offset  is  the 
difference  between  the  image  time  mark  and  a  straight  line  fit  to  the  time  marks. 

3.4  Modulation  Transfer  Function  of  the  System 


The  modulation  transfer  function  (NiTF)  of  an  optical  system  describes  the  effect  of 
aberrations  and  diffraction  on  an  image  formed  through  the  system.  If  /,(£)  is  the 
image  intensity  distribution  expected  from  aberration-free  geometrical  optics,  then 
the  real  image  intensity  distribution  /r(F)  is  given  by, 

lAD  =  .F-‘{MTF(it)  •  .r{/,(£)}}  (3.1) 

The  attenuation  of  the  spectral  power  is  given  by  the  square  of  the  MTF 

functions. 

The  components  which  contributed  significantly  to  the  total  MTF  are  the  entrance 
pupil,  the  detector  pixel  shape,  and  the  earth’s  atmosphere.  The  total  MTF  is  the 
product  of  the  MTF's  of  these  individual  components.  The  atmospheric  MTF  will  be 
treated  in  section  5.4.1. 
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3.4.1  Entrance  Pupil 

The  entrance  pupil  is  defitiefl  b}’  the  diameter  d  (48  cm)  of  the  heat  shield  diaphragm 
in  front  of  the  objective.  The  MTF  of  the  pupil  is  given  by  the  autocorrelation  of  the 
pupil  function  (a  disk).  For  wavenumbers,  0  <  k  <  t*Nyq, 

.\ITF(A-)  =  -  [cos-‘  K  -  /cn/T^I  (3.2) 

^  (  if  pixels  k 
w  leie  K  ^  field  of  view  /  A'Nyq 

3.4.2  CCD  Camera 

The  effect  of  summing  light  over  the  area  of  a  CCD  pixel  is  equivalent  to  convolving 
the  image  distribution  with  the  pi.vel  shape,  and  sampling  the  new  distribution  at 
points  centered  on  each  pixel.  Thus,  the  MTF  of  the  detector  is  the  modulus  of  the 
Fourier  transform  of  the  pixel  function  (a  square). 

MTF(A)  =  sincKx  sincKy  (3.3) 

where  k  =  --j - 

2  A'Nyq 

sin  T.r 

and  smc,r  = - 

■KX 

This  can  be  approximated  by  a  cylindrically  symmetric  function  with  no  more  than 
a  2%  error. 

MTF(A)  =  sincK  (3.4) 

The  square  of  the  MTF  functions  are  plotted  in  figure  3.12  for  the  entrance  pupil  and 
CCD  detector. 
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Chapter  4 


Observations 


The  original  topic  of  this  dissertation  was  "sunspot  seismology,”  that  is,  the  study 
waves  propagating  into  and  out  of  sunspots.  That  topic  was  set  aiside  when  waves  were 
discovered  at  siiljstantially  higher  wavenumber  and  frequency  than  had  previously 
been  observed.  The  new  topic  of  analysis  became  the  characterization  of  properties 
of  these  waves.  Many  of  t  he  observing  runs  are,  therefore,  extraneous  to  the  analysis 
described  here.  They  have  Ijeen  included  in  the  list  of  data  in  this  chapter,  but  no 
calibration  or  reduction  has  been  performed  on  them. 


4.1  Observing  Procedure 

The  observing  runs  consist  of  a  repeating  sequence  of  filtergrams  with  a  temporal 
spacing  of  about  3.5  s  between  filtergrams.  The  wavelength,  filter  width,  polariza¬ 
tion,  and  relative  exposure  time  for  each  filtergram  in  the  sequence  is  defined  in  the 
sequence  file  residing  on  the  C'OMP.AQ  computer  that  runs  the  filter. 

Calibration  data  was  taken  before  and  after  observing  runs.  The  order  of  calibra¬ 
tion  and  observations  was  as  follows: 

1.  Dark  Currents:  4  measurements  of  dark  currents  for  each  exposure  setting. 
Done  at  the  beginning  and  end  of  each  day. 
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2.  Line  Centering:  Definition  of  the  nominal  line  center  from  which  wave¬ 
length  offsets  are  measured.  The  photometer  detector  was  used  to  inte¬ 
grate  light  over  the  field  of  view  with  the  target  region  centered  in  the 
field. 

3.  Slide  Projector  Flats:  4  repetitions  of  flat  field  mecisurements  using  slide 
projector  some  e. 

4.  Observing  Run:  One  or  more  tapes-worth  of  observations  of  target  region. 

5.  Solar  Flats:  4  lepetitions  of  flat  field  measurements  using  defocused  sun¬ 
light  source. 

6.  Slide  Projector  Flats:  Repeat  of  step  3. 

The  line  centering  step  is  performed  by  scanning  the  filter  through  a  range  of 
wavelengths  (from  short  to  long),  and  plotting  intensity  versus  wavelength.  This  step 
was  not  performed  correctly  due  to  a  v?sign  flaw  in  the  photometer  interface  with 
the  COMPAQ  computer.  The  computer  did  not  wait  for  the  full  signal  integration 
time  before  reading  the  light  intensity.  In  the  resulting  line  profile  scan,  the  intensity 
points  lagged  behind  the  wavelength  points.  This  caused  the  line  to  appear  shifted 
to  longer  wavelength.  Tlie  line  was.  therefore,  “centered”  to  an  erroneously  short 
wavelength.  This  offset  of  — 10  m.A  for  the  Na  Di  line  was  discovered  by  making  a 
line  scan  using  the  CCD  camera  instead  of  the  photometer  (see  figure  4.1). 


4.2  Calibration  Method 

The  purpose  of  calibration  is  to  iiroduce  images  with  a  signal  to  illumination  ratio 
that  is  independent  of  both  the  pi.xel  location  and  the  wavelength  setting  of  the 
tunable  filter.  The  procedure  described  here  is  more  complicated  than  the  standard 
one  used  by  the  Lockheed  group.  This  is  necessary  to  provide  an  accurate  calibration 
for  dopplergrams.  The  flat  fielding  procedure  calibrates  the  whole  system,  optics  and 
detector  together. 

The  CCD  is  a  very  linear  device,  provided  the  light  intensity  stays  below  saturation 
levels.  Response  of  a  pixel  can  be  characterized  by  a  gain  and  am  offset.  The  offset 
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Figure  4.1:  Line  scan  of  .Xa  Di  Line  using  the  CCD  camera.  Signal,  integrated  over 
the  field  of  view,  vs.  wavelength  offset. 

is  determined  by  the  amount  of  charge  that  accumulates  on  the  CCD  in  the  absence 
of  illumination  (the  “dark  current”).  It  must  be  measured  for  each  exposure  time. 
Each  pixel  has  a  characteristic  gain  and  dark  current.  The  tunable  filter  is  also 
responsible  for  some  variation  of  the  signal  strength  as  a  function  of  wavelength  A 
and  spatial  position  x.  The  net  response  of  the  system  is  to  produce  filtergrams  Fi{x) 
at  wavelengths  Aj,  given  a  source  intensity  distribution  S(x,  A). 

F.(x)  =  JS(x,X)Ai(x,\)d\  +  D(x)  (4.1) 

where  y4i(x,A)  is  the  product  of  filter  transmission  and  CCD  gain  at  filter  setting  A^, 
and  D(x)  is  the  CCD  dark  current. 

Because  the  filter  bandpass  is  narrow  compared  with  spectral  features  in  the 
source,  we  can  approximate: 

J S(x,  A)  >l,(f, A)  dX  =  S(f, A.)  J Ai(x,  X) dX.  (4.2) 


Thus 


F,(x)  =  S(f,AOB<(f)  +  D(x) 


(4.3) 
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where  S,(f)  =  jAi(x,X)dX.  (4.4) 

In  order  to  determine  .5'.  tlie  quantity  of  interest,  we  must  first  measure  5,  and  D. 
The  dark  current  D  is  easily  measured  by  making  filtergrams  with  no  illumination 
source.  For  5  =  0,  we  ha\e  D  =  Fi. 

If  we  had  an  illuminaiion  source  that  were  both  flat  in  spectrum  and  flat  across 
the  field  of  view,  we  could  measure  Bi  directly.  The  sources  available,  however,  were: 

1.  Defocused  Sunlight  Flat  across  field,  but  not  spectrally  flat. 

2.  Slide  Projector  Lamp  Sl(  v):  Flat  in  spectrum,  but  not  flat  across  field. 


The  technique  is  to  make  one  filtergram  Ff  at  wavelength  A,,  in  the  continuum 
(away  from  the  absorption  line),  using  defocused  sunlight.  And,  to  make  filtergrams 
F^  at  Ac  and  Ff'  at  A,  using  the  slide  projector  lamp.  This  gives  us  the  following 
three  measurements: 

Ffix)  =  55(Ac)  Bcix)  +  D5(x)  (4.5) 

Ff'ix)  =  .?£(.?)  Bc(.?)  +  Dl(x)  (4.6) 

F!'(x)  =  5t(.r)  Bi(£)  +  Dl(x)  (4.7) 

Given  the  dark  currents  Di  aiul  D5.  we  can  solve  for  the  gain  correction  Qi  at  Aj 
defined  as  ^.(.r)  =  fl,(,r)”‘. 


Q,(F) 


[Ft  -  Dl] 

[Ft  ~  DiWFf  -  Ds] 


(4.8) 


There  is  some  arbitrary  overall  normalization  which  depends  on  5s(Ac).  In  practice 
the  spatial  average  of  Q,  at  the  A,  filter  position  was  defined  as  unity,  and  the  other 
gain  corrections,  normalized  accor<lingly. 

With  the  gain  corrections  and  dark  currents  in  hand,  we  can  determine  arbitrary 
source  distributions  from  rilt«’rgrams. 


S{x.X.)  =  g.(F)[F.(x)-D(x)] 


(4.9) 
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4.3  Solar  Observations 

The  sequence  of  filtergrains  for  the  sunspot  seismology  study  was  named  SEISNAFE. 
It  was  designed  to  determine  the  following  quantities: 

1.  Doppler  shift  ir  mid-chromosphere  (in  sunspot  umbrae). 

[H  I  (Balmer)  o — line  core] 

2.  Doppler  shift  in  photosphere. 

[Na  Di — wings] 

3.  Doppler  shift  in  low  chromosphere. 

[Na  Di — core] 

4.  Doppler  shift  in  photosphere. 

[Fe  I  (557(i.-\)] 

5.  Continuum  intensity  in  green  and  red. 

[55 75. A.  5StM.\! 

6.  Photospheric  magnetic  field  (line  of  sight  component). 

[Fe  I  (().302.\)) 

The  sequence  repeated  e\er\'  17  s.  Filter  parameters  appear  in  table  4.1.  The  corre¬ 
sponding  calibration  sequences  were  .N'.AFEFL.AT  for  solar  flats,  and  SEISSP  for  slide 
projector  flats.  The  observing  runs  of  appreciable  length  are  listed  in  table  4.2. 

The  sequence  of  filtergrains  used  for  the  high-frequency  oscillation  study  was 
named  NAF.AST.  It  was  designed  to  make  dopplergrams  in  the  low  chromosphere 
and  photophere.  using  filtergrains  in  the  core  and  wings  of  Na  Di,  respectively.  All 
other  measurements  wer<-  eliminated  to  increase  the  repetition  rate  to  once  every 
17.7  s.  The  sequence  parameters  appear  in  table  4.3  and  the  corresponding  calibra¬ 
tion  sequences  were  .\.A1-  F.F1..\T  for  solar  flats,  and  N.AFASTSP  for  slide  projector 
flats. 

The  NAF.AST  sequence  was  run  only  once  on  .August  2.  1990  beginning  13:51 
UT.  The  target  region  was  an  active  region  (.AR61S1)  containing  a  small,  decaying 
sunspot.  This  is  shown  on  a  fuii-disk  magnetogram  in  figure  4.2.  The  run  lasted 
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Table  4.1: 

:  Filter 

!^arameters  for  SEISNAFE  Sequence. 

Line  (.4) 

(ill. A) 

Filter  Mode 
(width) 

Polarization 
(0  =  linear) 

Purpose 

H  I  a  6562.8 

-400 

0 

400 

wide 

0 

dopplergram 

Na  D,  .5895.9 

-10 

0 

40 

narrow 

0 

dopplergram 

Na  Di  .5895.9 

-175 

-125 

125 

175 

narrow 

0 

dopplergram 

Na  Dj  5895.9 

-2000 

narrow 

0 

continuum 

Fe  1  .5576.1 

-00 

-:10 

0 

40 

90 

wide 

0 

dopplergram 

Fe  I  .5576.1 

-750 

wide 

0 

continuum 

mmimgii 

-60 

-60 

narrow 

LCP 

RCP 

magnetogram 
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Table  4.2:  La  Palma  Observing  Runs  using  SEISNAFE  Sequence. 


Date 

(1990) 

Time 

(UT) 

Duration 

(min) 

Seeing 

FOV 

Target  Region 

28  Jul 

my 

■n 

goocl/v.good 

75" 

AR6168  (medium  sunspot) 

Bai 

fair/var. 

150" 

AR6168 

29  Jul 

156 

fair/var. 

150" 

AR6174  (medium  sunspot) 

126 

goocl/var. 

150" 

30  Jul 

7:50 

281 

variable 

224" 

13:12 

293 

variable 

224" 

31  Jul 

8:34 

237 

gootl/var. 

224" 

—  flare  in  morning 

13:09 

347 

variable 

224" 

1  Aug 

8:08 

256 

goocl/var. 

224" 

AR6181  (small  sunspot) 

ISQ 

13:35 

3:o 

goocl/fair 

224" 

AR6179  (big  sunspot) 

HSil 

7:46 

151 

goocl/fair 

150" 

—  overexp’d  for  umbral  osc. 

2  Aug 

10:52 

101 

lair/ var. 

224" 

AR6186  (big  sunspot  on  limb) 

5  Aug 

9:20 

100 

fair/ vav. 

224" 

c\uiet  region  near  .\R6192 

Table  1.;: 

;  Filter 

'Parameters  for  N.AF.AST  Sec 

uence. 

Line  (A) 

Olfset 

(niA) 

Filter  .Mode 
(width) 

Polarization 
(0  =  linear) 

Purpose 

Na  Di  5895.9 

-10 

0 

10 

narrow 

0 

dopplergram 

Na  Di  5895.9 

-175 

- 1 25 

1 25 

175 

narrow 

0 

dopplergram 

.Va  D|  5895.9 

-2000 

narrow 

0 

continuum 
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Figure  4.2;  Full-disk  iua«iietv)j>,iain  from  the  National  Solar  Observatory  at  Kitt  Peak. 
The  field  of  view  containing  .\H6181  is  outlined. 

143  min  with  good/varialile  seeing  that  deteriorated  in  the  last  30  min  of  the  run.  To 
obtain  more  information  about  the  active  region,  the  SEISN.AFE  sequence  was  run 
before  and  after  the  .\.\F.\S  1  runs.  .A  continuum  image  is  shown  in  figure  4.3,  a  high 
resolution  magiu’logram.  in  lignie  1.1.  and.  an  offband  H  q  image,  in  figure  4.5. 
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Figure  4.3:  Continuum  intensity  (5575A)  image  of  AR6181.  The  small  tick  marks 


Figure  4.4:  Magnetogram  of  .AR61S1.  Network  field  is  visible  around  the  sunspot. 
The  small  tick  marks' are  5". 


CHAPTER  4.  OBSERVATIOyS 


60 


Figure  4.5:  H  a  image  of  .AR61S1  (400  niA  in  the  blue  wing).  The  small  tick  marks 
are  5". 


Chapter  5 

Data  Reduction 


The  raw  data  used  for  this  dissertation  reside  on  a  (full)  2.2  Gigabyte  Exabyte  tape.  In 
producing  a  power  spectiuni.  the  volume  of  data  is  reduced  by  a  factor  of  ~  3  x  10"*. 
There  is  more  to  the  reduction  procedure,  however,  than  extracting  a  manageable 
amount  of  information  from  a  cumbersome  data  set.  The  quality  of  the  data  depends 
on  careful  calibration  to  account  for  instrument  response,  image  rotation  and  drift 
caused  by  telescope  motion,  and.  distortion  caused  by  the  atmosphere.  The  absolute 
spatial  scale  of  the  images  must  also  be  determined  accurately.  This  chapter  describes 
the  details  of  and  cross-cliecks  for  the  data  reduction  procedures. 

The  reduction  steps  are  listed  below  in  the  order  they  were  performed: 

1.  Flat  Fielding:  Calibrate  for  variation  in  filter  transmission  and  CCD  gain, 
and  subtract  CCD  dark  current. 

2.  CCD  Flaw  Correction:  Interpolate  over  nonfunctional  pixels. 

3.  Correction  for  Atmospheric  Blurring:  Normalization  of  measured  MTF. 

4.  Rigid  Rotation:  Compensate  for  image  rotation  caused  by  telescope  track¬ 
ing  motions  during  observations. 

5.  Correction  for  Atmospheric  ".Shimmer.”  (Destretching) 

6.  Calculation  of  Dopplergrams  and  Line  Center  Intensity. 

7.  Rigid  .Alignment:  Compensate  for  image  drift  caused  by  telescope  tracking 
errors. 


(>1 
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8.  Fourier  Tiaiisf'i)im  of  Dopplergrams  and  Intensities:  Described  in 
section  2.4. 

9.  .-\ngiilar  Scale  Determination:  Convert  units  fronr  pixels  to  arcseconds. 

10.  Distance  Scale  Determination:  Measurement  of  foreshortening  and  direc¬ 
tion  toward  disk  c<’nter.  Described  in  section  6.1.1. 

11.  Measurement  of  Oscillation  Properties:  Described  in  chapters  6  and  7. 


5.1  Flat  Fielding 


The  flat  fielding  scheme  is  explaitietl  in  section  4.2.  W’e  can  use  equation  4.9  to 
calibrate  filtergrams.  given  dark  currents  D  and  gain  corrections  Qi.  The  dark  current 
is  an  average  of  four  filtergrams  made  with  a  light  stop  placed  in  front  of  the  filter. 
The  gain  correction  is  coiistnufed  using  equation  4.8  as  follows. 

The  three  flat  field  measurements  Ff .  and  Ff'  are  each  made  from  an  average 

of  four  filtergrams.  Because  the  filter  stepping  motors  did  not  always  find  the  correct 
positions,  it  was  sometimes  neces.saiy  to  throw  out  one  of  the  four  measurements 
and  use  only  three  to  make  the  average.  .A  filtergram  was  thrown  out  if  the  spatial 
average  of  the  central  64  x  61  pi.xels  of  the  image  differed  by  more  than  5<t  from  the 
spatial  averages  of  subareas  of  the  other  three  images,  where  a  is  defined  by  inefficient 
statistics. 


/, 


a,  = 


max 


mm 


.V, 


(5.1) 


The  gain  correction  is  the  product  of  (Ff  —  Ds)~'  and  (F/'  —  DL)f(Ff'  —  Dl). 
The  latter  is  the  ratio  of  the  transmittance  in  continuum  to  transmittance  an  A^. 
This  ratio  should  be  slowly  varying  over  the  field  of  view.  The  dark  currents  are 
not  constant  with  time,  however,  and  the  very  long  exposure  times  required  for  slide 
projector  flats  (2-3  s)  caused  accumulation  of  large  charges  at  some  small  groups  of 
pixels.  It  was,  thus,  necessary  to  "deflaw"  F/”  and  Ff"  (see  next  section),  and  smooth 
the  transmittance  ratio  using  a  Gaussian  point-spread  function  of  18  pixel  fwhm. 

The  procedure  described  here  produced  (flat  fielded)  images  with  signal  to  illu¬ 
mination  ratios  that  wen*  constant  to  within  about  3%  across  the  field  of  view,  and, 
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more  importantly,  to  williin  aliout  o%  across  the  wavelength  range  of  an  absorption 
line.  These  numbers  \vei'<‘  measured  from  filtergrams  of  defocused  sunlight  (the  flat 
field  illumination  source)  at  continuum  wavelengths. 

5.2  CCD  Flaw  Correction 

Of  the  2.5  X  10®  pixels  (in  2  x  2  summing  mode),  4513  were  not  functional.  .An 
additional  491  pi.xels  had  vceptionally  high  dark  currents  for  exposure  times  greater 
than  about  500  ms.  To  hiinimize  prolilems  caused  by  the  missing  data,  values  were 
filled  in  at  these  pixel  locations  l)y  linear  interpolation  from  adjacent  pixels.  This 
procedure  is  termcxl  '■deflawing". 

A  good  map  of  nonfimci ional  pi.xel  locations  was  available  for  exposure  times 
below  500  ms.  For  the  longer  exposure  times,  a  map  was  constructed  using  the  dark 
current  runs.  To  do  this,  a  cubic  trend  was  fitted  to  each  row  of  a  filtergram.  Pixels 
which  differed  from  the  tieml  by  move  than  were  flagged  as  “bad”,  and  their 
locations,  entered  in  the  map. 


5.3  Image  Rotation 

The  altitude-azimuth  axes  of  mirrors  Ml  and  M2  cause  the  solar  image  to  rotate 
during  the  day  as  the  telescop**  tracks  the  sun.  The  rotational  direction  and  rate  are 
functions  of  <late  and  time  of  <lay.  and  can  be  calculated  from  geometrical  consider¬ 
ations.  Standard  software  exists  to  determine  the  angle  of  rotation  and  regrid  the 
images  by  bicubic  interpolation.  This  removes  the  apparent  rotation  in  a  sequence  of 
images,  but  some  data  is  lost  from  the  corners  of  the  images. 


5.4  Atmospheric  Distortion 

Turbulence  in  the  earth's  atmospher**  will  modify  wavefronts  of  light  before  they  reach 
the  telescope.  Any  plane  wa\e  passing  through  the  atmosphere  will  be  distorted. 
The  resulting  wavefront.  how<*ver.  can  be  approximated  by  a  set  of  planar  patches. 
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The  isoplanatic  angle  describes  the  angular  extent  over  which  image  distortions  are 
correlated.  For  excellent  seeing  conditions,  it  can  be  over  100"  (von  der  Luhe  1988). 
For  poor  seeing,  it  can  b<-  less  than  1". 

A  change  in  angle  of  incidence  of  a  planar  wavefront  produces  an  image  displace¬ 
ment.  When  different  parts  of  an  image  are  displaced  in  different  directions,  the 
image  appears  distorted.  The  distortions  change  with  the  atmospheric  turbulence. 
This  can  occur  over  a  range  of  timescales,  from  1  to  1000  ms,  with  most  power  in  the 
10-100  ms  range. 

5.4.1  MTF  Normalization 

The  distortion  of  a  wa\<‘front  from  a  given  object  point  differs  over  the  telescope 
aperture.  Moreover,  local  displacements  of  up  to  an  arcsecond  cause  the  image  to 
move  around  in  times  of  order  10  ms.  For  exposures  lasting  tens  to  hundreds  of 
milliseconds,  these  changing  distortions  will  be  superimposed  to  produce  a  blurred 
filtergram.  The  effect  is  ecpii valent  to  convolving  the  image  with  some  unknown 
point  spread  function.  The  Fourier  transform  of  the  point  spread  function  is  called 
the  MTF.  It  is  assumed  here  that  the  point  spread  function  is  the  same  everywhere 
in  the  field  of  view  of  a  given  liltergrain. 

The  atmospheric  MTl-  can  l)e  estimated  from  the  spatial  frequency  spectrum  of 
a  filtergram  S{f).  contour  plot  of  spatial  frequency  spectra  is  shown  in  figure  5.1 
for  the  —40  m.A  filtergram  set.  The  sharpest  filtergrams  in  the  set  are  assumed  to 
be  diffraction  limited  (atmospheric  MTF  =1).  Since  the  structure  of  solar  images 
should  not  change  apprecial)ly  over  the  observing  run,  the  spatial  spectra  of  diffraction 
limited  images  should  be  a  constant  >,«««(/)•  Any  difference  between  the  spectrum 
of  a  given  filtergram  and  must,  therefore,  be  due  to  the  atmospheric  MTF. 

.s'(/)  =  .MTF..„(/).5n,„(/)  (5.2) 

In  practice  the  log-spectra  were  compared  with  the  mean  log-spectrum  log 5. 
Log-spectra  of  individual  liltergrams  were  smoothed  by  convolution  with  a  cosine-bell, 
and  the  difference  betwer’ii  log and  the  smoothed  spectrum  was  used  as  a  correction 


Lo9(RMS  amplitude) 
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Figure  5.1:  Spatial  frequency  log-spectra  of  Na  Di  —40  mA  filtergr2uns.  The  solid 
central  line  is  the  mean  log-power:  the  dashed  lines  are  la  contours;  the  outer  dotted 
lines  are  the  maximum  and  minimum  power  at  a  given  frequency  of  any  image. 
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factor  to  the  filtergrani  freciuencies.  Thus,  the  "MTF  normalization’"  procedure  forced 
the  spatial  spectrum  of  each  filtergrani  to  look  like  the  average  spectrum. 

In  principle  the  individual  filtergrams  could  have  been  normalized  to  5max-  This 
would  have  involvetl  very  large  correction  lactors  for  the  poorer  filtergrams,  and  would 
have  boosted  the  noise  in  t  iiese  filtergrams  considerably.  Noise  in  filtergrams  produces 
noise  in  dopplergrams  aiul  line  center  intensities  in  a  nonlinear  manner.  So,  instead, 
the  difference  between  the  mean  spectrum  and  was  retained  as  a  correction 

factor  to  be  applied  to  the  oscillation  power.  (See  chapter  6.)  As  a  final  caveat,  it 
should  be  noted  that  even  the  best  images,  assumed  above  to  be  diffraction  limited, 
actually  include  some  atmospheric  blurring.  Hence,  the  correction  factor  will  provide 
an  underestimate  lor  po\V(>r  at  the  highest  wavenumbers. 

The  MTF  tiortiialization  procedure  was  checked  by  comparing  filtergrams  taken 
in  poor  seeing  with  filtergrams  taken  at  almost  the  same  time  in  good  seeing.  .After 
normalization,  both  filteigratns  showed  the  same  solar  structures. 

5.4.2  Destretching 

The  displacemetits  cause<l  b\  atmospheric  turbulence  vary  over  the  field  of  view. 
They  are  independent  of  one  another  over  angular  scales  greater  than  the  isoplanatic 
angle.  Relative  displacentetits  between  successive  images  can  be  measured  by  local 
correlation  tracking  (.\o\ (‘tuber  1986).  This  procedure  finds  the  displacement  of  a 
small  subarea  of  one  image  with  respect  to  the  other  from  the  cross-correlation  of  the 
two.  A  displacement  vector  field  is  measured  using  a  patchwork  of  subareas  of  the 
images.  .\  relative  disphn  ement  grid  of  64  x  64  vectors  was  made,  correlating  each 
filtergram  with  the  preceding  one.  This  was  done  using  standard  software. 

Local  displacements  are  present  not  only  because  of  seeing,  but  also  because  of 
wave  motion  and  horizontal  flows  on  the  sun.  VVe  wish  to  remove  the  seeing  displace¬ 
ments,  but  keep  the  real  solar  motions.  It  is  po.ssible  to  distinguish  the  two  on  the 
basis  of  timescale  of  the  motion.  MorizoiUal  flows  have  characteristic  timescales  of 
hours,  oscillations  have  timescales  of  minutes,  and  seeing  has  a  timescale  of  tens  of 
milliseconds.  For  filtergrams  taken  seconds  apart,  seeing  motions  will  appear  random, 
while  solar  ntotions  will  be  coherent. 
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The  seeing  displacement  grid  was  determined  from  the  relative  displacement  grid 
as  follows.  For  each  vector  component  at  each  grid  location,  the  relative  displacements 
of  preceding  filtergrams  were  summed  to  produce  a  total  net  displacement  (from  the 
first  frame).  The  total  displacement  was  smoothed  by  convolution  with  a  Gaussian 
of  90  s  (5  bin)  fwhm,  and  this  smooth  trend,  subtracted  to  produce  the  seeing  grid. 

It  is  important  to  hav<-  high  temporal  resolution  in  the  data.  Correlation  tracking 
can  lock  onto  the  crests  of  waves.  So.  if  the  impulse  response  function,  used  to  smooth 
the  data,  is  wider  than  ~  120  s,  the  oscillations  will  be  removed  from  the  data  along 
with  the  seeing.  The  standard  l.ockhewl  method  for  finding  the  seeing  grid  is  designed 
for  horizontal  flow  studies.  It  uses  boxcar  averages  of  20  or  more  bins.  This  would 
be  inappropriate  for  seismology  studies. 

With  the  seeing  displacements  measured,  the  images  can  be  “destretched”  to 
their  undistorted  states.  I'he  algorithm  to  do  this  treats  images  as  rubber  sheets,  and 
stretches  the  sheets  (by  regridding)  based  on  the  local  displacement  vectors.  This  part 
of  the  process  also  used  standard  software.  The  set  of  filtergrams  for  each  wavelength 
was  destretched  independently  of  the  other  wavelength  sets.  Destretched  filtergrams 
taken  at  the  same  time,  but  <Iifl'erent  wavelengths,  were  then  visually  compared  to 
ensure  their  features  lined  up. 


5.5  Dopplergrams  and  Intensities 

Velocity  in  m/s  and  line  center  intensity  in  arbitrary  units  were  calculated  at  each 
pixel  using  equations  2.2  and  2.d.  Because  of  the  line  centering  problem  described  in 
section  4.1,  the  line  center  velocity  and  intensity  were  calculated  using  the  —125,  —40, 
and  +40  inA  filtergrams.  The  actual  offsets  of  these  measurements  from  the  average 
line  center  was  about  —<'>5.0.  and  +80  m.4.  Each  dopplergram  and  intensity  image 
was  given  a  time  mark  a\eraged  from  the  time  marks  of  the  three  filtergrams. 

Filter  tuning  errors  and  seeing,  which  deteriorated  toward  the  end  of  the  run, 
produced  some  ijoor  quality  filtergrams  in  the  data  set.  This  in  turned  caused  some 
nonsensical  velocity  and  intensity  determinations.  Pixel  velocities  that  were  greater 
in  magnitude  than  10  km/s  and  intensities  that  were  a  factor  of  10  greater  than 
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average  were  replaced  with  tlie  mean  velocity/iiitensity  for  that  frame. 


5.6  Image  Drift 

The  rotational  regriddiiig  assumes  the  images  were  rotating  about  a  central  point. 
In  fact  they  rotate  about  the  [joint  centered  in  the  spot  tracker.  This  inconsistency 
produces  a  translational  motion  of  the  solar  image.  The  point  (on  the  sun)  centered 
in  the  spot  tracker,  moieo\  er,  changes  with  time.  This  is  caused  by  variable  weighting 
of  parts  of  the  active  region  as  the  region  evolves  and  as  seeing  conditions  change. 

From  equation  ‘2.27.  it  is  evident  that  to  keep  oscillation  frequencies  from  shifting 
by  more  than  60  //Hz  (0.')  lr<  (|iieiK  V  bins)  at  a  wavenumber  of  3  .\Im~'.  the  image 
drift  must  be  less  than  120  m/s.  This  is  equivalent  to  1  !  pixel/hr.  Image  drift  is 
measured  from  the  cross-/ orrelat  ion  of  successive  images.  (It  is  defined  as  the  distance 
of  the  maximum  of  the  cross-correlation  from  the  origin.) 

The  image  drift  is  trarefl  out  o\er  fijjie  in  figure  5.2.  It  was  measured  using  the 
central  64  x  61  pi.xels  of  every  t/uith  —125  m.\  filtergrnm.  .As  can  be  seen  from  the 
figure,  drift  can  be  measured  to  significantly  better  than  1  pixel.  The  meaisured  drift 
w'as  smoothed  by  cotivoh  itig  I  Ik*  lime  seri<^  of  each  component  of  the  vector  with  a 
Gaussian  of  6  l/in  Ju'lnn.  I  lie  tread  was  then  interpolated  to  the  dopplergram  ®nd 
intensity  time  marks  aiul  renuned  from  these  images  by  regridding. 


5.7  Angular  Scale  Determination 

For  full  disk  images  of  the  sun.  the  image  scale  can  be  determined  by  measuring 
the  radius  of  the  solar  disk,  l  or  high  r/’soiution  images,  however,  the  full  disk  is 
not  visible,  and  indir/’ct  ineth/uls  must  be  used.  The  scale  of  the  La  Palma  images 
(0.8x  lens)  was  determined  b\'  comparison  of  magnetograms  from  La  Palma  and  the 
Kitt  Peak  National  Observatory  (figures  4.4  and  4.2).  The  Kitt  Peak  magnetogram 
was  made  at  15:17  I’T  on  .August  2.  1!)!)0.  and  has  a  scale  of  1"/ pixel  to  within  1% 
(Harvey  1901).  The  La  Palma  mngn/’togram  was  interpolated  from  magnetograms 
made  at  13:35  FT  and  l(i;27  FT. 
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Figure  5.2;  Image  drift  during  observing  run.  Image  motion  determined  by  cross¬ 
correlation  is  traced  out  for  the  observing  run.  Points  represent  iiitergrams  taken  177 
s  apart;  solid  line  shows  low-pass  filter  of  points;  circles  represent  18  min  time  marks 
(<  =  0  is  at  upper  right). 
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The  interpolation  iiietliod  accounted  for  evolution  of  the  active  region  by  measur¬ 
ing  the  motion  of  the  magnetic  field  structures  using  local  correlation  tracking.  The 
two  magnetograms  were  t  In-n  stretched  by  the  time-weighted  local  displacement  grid. 
Finally,  the  time-weighted  signals  of  the  stretched  magnetograms  were  added. 

To  compare  magnetogi  ams.  -13  points  on  the  La  Palma  magnetogram  were  flagged 
manually  by  cursor  on  a  P1.\.\R  cou\puter.  Corresponding  points  were  flagged  on  the 
Kitt  Peak  magnetogram.  Distances  between  a  point  and  the  centroid  of  the  43  points 
were  then  compared  for  t  he  two  magnetograms.  The  field  of  view  for  the  La  Palma 
images  was.  thus,  determined  to  be  {'i'i-L.jy  ±  0.70)  Kitt  Peak  pixels.  .Assuming  the 
Kitt  Peak  scale  measurement. 

La  Palma  Field  of  View  =  224"  ±2"  (square)  (5.3) 


f* 


Chapter  6 

Measurements  from  Power 
Spectra 


The  sets  of  dopplergiams  and  intensities  were  apodized  in  the  space  and  time  di¬ 
rections  with  Blackman- Harris  functions  (equation  2.32).  Three-dimensional  FFT’s 
were  calculated  for  the  full  field  of  view,  as  well  as  for  the  three  79"  x  79"  subfields 
shown  in  figure  6.1.  The  temporal  and  spatial  Nyquist  frequencies  were  28.3  mHz  and 
9.9  Mm“‘,  respectively.  The  frequency  bin  size  was  117  ^Hz,  and  the  wavenumber 
bin  size,  0.039  Mm~'  for  the  full  field,  and,  0.11  Mm"*  for  the  subfields. 

The  squared  modulus  of  each  Fourier  transform  was  calculated  to  produce  a  3-d 
power  spectrum  P(k,i/).  The  dispersion  relations  given  by  equations  1.7  and  1.9 
are  symmetric  about  the  A-  =  0  axis.  They  define  “trumpet”  shaped  surfaces  (Hill 
1988)  of  approximate  form  //  oc  jA'I'/*.  (See  figure  6.2.)  This  synunetry  allows  us 
to  integrate  the  power  o\er  the  azimuthal  coordinate,  in  order  to  increase  the  signal 
to  noise  ratio.  The  resultant,  collapsed  power  spectrum  P{k,i/)  will  show  ridges  of 
power  at  frequencies  given  by  the  dispersion  relations  for  the  p-modes  auid  /-modes. 

Properties  of  the  oscillations,  such  as  the  dispersion  relation  and  mode  power,  were 
measured  from  the  collapsed  power  spectra.  This  was  accomplished  by  fitting  one- 
dimensional  slices  of  the  s|)ectra.  Fits  were  then  compared  with  theoretical  predictions 
to  evaluate  models  of  the  solar  atmosphere.  Magnetic  versus  nonmagnetic  r^ions 
were  compared  to  search  for  rlifferences  in  wave  properties  between  the  two  regions. 
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Figure  6.1:  Magnetogram  of  .\R61S1.  The  tlu’ee  outlined  subfields  were  used  sepa¬ 
rately  to  measure  wave  properties  in  network  field  and  quiet  regions.  The  small  tick 
marks  are  -5". 
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6.1  Collapsing  3-d  Power  Spectra 

The  azimuthal  symmetry  of  equations  1.7  and  1.9  is  realized  only  in  solar  coordinates, 
comoving  with  the  gas.  In  practice,  one  or  both  of  these  conditions  may  not  be  met. 
The  frequency  shift  described  by  equation  2.27  occurs,  with  a  cosine  dependence  on 
azimuth  0,  if  there  is  a  rela  tive  (horizontal)  motion  between  the  observer’s  coordinates 
and  the  gas.  The  rigid  alignment  procedure  described  in  section  5.6  keeps  the  observa¬ 
tional  coordinates  fixed  with  respect  to  the  sunspot.  It  has  been  shown,  however,  that 
there  is  relative  motion  bc'tween  magnetic  flux  and  nonmagnetized  gas  (Snodgrass  & 
Ulrich  1990).  Hill  (19SS)  finds  flows  of  ~  100  m/s  below  the  photosphere. 

Visual  inspection  of  the  power  spectra  did  not  reveal  evidence  of  flows.  Thus, 
if  present,  flows  produce  frequency  shifts  of  less  than  one  frequency  bin  (117  /iHz). 
If  plane  waves  propagated  isotropically,  a  sinusoidal  frequency  shift  with  0  would 
serve  to  widen  the  ridges  in  the  collapsed  power  spectra.  In  the  measured  spectra, 
however,  power  depends  mildly  on  angle  0.  Thus,  a  horizontal  flow  would  produce  a 
slight  asymmetry  in  addition  to  broadening  a  ridge. 

The  second  condition  affecting  the  symmetry  of  the  trumpets  is  geometrical.  For 
observations  away  from  disk-center,  the  projection  of  solar  coordinates  onto  obser¬ 
vational  coordinates  caus«'S  foreshortening.  Waves  propagating  toward  or  away  from 
disk-center  appear  with  shorter  wavelength  (higher  wavenumber)  than  waves  prop¬ 
agating  in  the  perpendicular  direction.  This  foreshortening  distorts  the  cylindrical 
symmetry  of  the  trumpets  to  give  them  an  elliptical  rather  than  circular  cross  section. 
Procedures  to  measure  and  compensate  for  this  effect  are  described  below. 

6.1.1  Distance  Scale  Determination 

The  distance  scale  of  the  images  or.  equivalently,  the  wavenumber  scale  of  the  power 
spectra  follows  from  trigonometry,  once  the  magnitude  and  direction  of  foreshortening 
have  been  determined.  This  measurement  was  made  by  fitting  an  ellipse  to  the  cross 
section  of  a  trumpet  in  a  3-d  power  spectrum. 

A  constant  frequency  slice,  in  tlie  5-minute  band,  of  the  full-field  Doppler  spectrum 
is  shown  in  figure  6.3.  The  slice  has  been  remapped  to  cylindrical  polar  ooturdinates. 
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Figure  6.3;  Constant  fiec|iiency  slice  of  the  3-d  power  spectrum.  A  slice  at  3.99  mHz 
has  been  remapped  to  c\  liiKlrical  polar  coordinates.  Log(poirer)  is  shown  as  a  function 
of  wavenumber  and  azimuth. 

So,  the  cross  section  of  a  trui\i|)et  appears  at  roughly  constant  wavenumber  for  all 
values  of  0.  These  horizontal  l)aiKls  of  power  are  the  “rings”  described  by  Hill  (1988). 
In  polar  coordinates,  an  e  llipse  centered  at  the  origin  has  the  form, 

k  =  Ao[l  -  (1  -/«^)  cos*(0  -0o)]~5  (6.1) 

where  ko  is  half  the  minor  axis  (the  true  wavenumber):  is  the  angle  of  the  major 
axis  (direction  toward  disk-center):  and  //  is  the  ratio  of  minor  to  major  axes  (the 
fractional  foreshortening).  close  inspection  of  figure  6.3  shows  that  each  of  the 
bands  exhibits  this  oscillatorx  behavior  of  k  as  a  function  of  $. 

The  values  of  ko,0u.  and  //  were  fitted  to  the  measured  wavenumber  k{0)  for  the 
Pi  ring  in  the  frequency  slice  //  =  3.99  mHz.  The  function  k{9)  was  defined  as  the 
wavenumber  at  peak  power  in  the  constant  frequency  slice.  It  was  determined  by 
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quadratic  interpolation.  I'he  measured  fractional  foreshortening  was, 

/t  =  0-9.56 

This  corresponds  to  a  distance  from  disk-center  of  17®,  which  agrees  with  an  approx¬ 
imate  mea-surement  made  from  the  full-disk  magnetogram  in  figure  4.2. 

The  assumption  of  a  planar  atmosphere  hais  been  made  here,  in  treating  the  fore¬ 
shortening  as  constant  across  the  field  of  view.  If  the  foreshortening  varies  appreciably, 
modes  of  different  wavenumber  will  be  mixed  in  the  3-d  power  spectrum.  This  will 
smear  the  trumpet  surfaces  in  {k\u)  space  in  a  manner  which  cannot  be  corrected. 
For  large  fields  of  view,  the  images  must  be  mapped  to  “flattened”  coordinates  before 
transforming  to  Fourier  space.  The  aplanarity  for  the  data  analyzed  here  was  not 
great  enough  to  warrant  such  treatment. 

The  224"  field  of  view  represents  13®  along  a  great  circle  on  the  Sun.  Since  the 
average  distance  from  disk-center  was  about  17®,  plane  waves  propagating  along  the 
direction  of  foreshortening  will  l)e  mixed  over  a  wavenumber  range  of  ±3.5%.  While 
this  may  seem  large,  only  the  waves  propagating  exactly  along  a  diameter  of  the  disk 
experience  mixing  to  this  extent.  Waves  propagating  in  the  perpendicular  direction 
are  not  mi.xed  at  all.  The  net  widening  of  the  ridges  is,  thus,  small.  This  can  be  seen 
in  figure  6.3.  The  width  of  the  ridges  api>ears  independent  of  0.  The  width,  therefore, 
is  not  an  artifact  of  the  coordinate  system. 


6.1.2  Azimuthal  Projection 

The  power  was  integrated  around  ellipses  in  the  (kx,ky)  plane.  To  do  this,  ko  was 
computed  from  equation  6.1  for  each  ordered  pair  (kx,kp)  in  the  3-d  spectrum,  and 
power  from  that  bin  was  added  to  the  appropriate  k  bin  in  the  2-d  spectrum.  Bins 
in  the  2-d  spectrum  were  norntalized  by  (multiplying  by)  the  ratio  of  the  area  of  an 
elliptical  band  to  the  number  of  bins  accumulated  from  the  3-d  spectrum.  The  area 
A  of  an  elliptical  band  was  a|}|>roximated  by, 

.4  =  «•— \/a*  -h  IP 


(6.2) 
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Figure  6.4:  Scaled  velocity- power  spectrum.  A  portion  of  the  collapsed  power  spec¬ 
trum  for  the  full  field  of  view,  with  variable  gray-scaling  in  the  k  and  t'  directions. 
Background  has  been  subtracted.  The  noise  visible  at  the  top  of  the  spectrum  is  an 
artifact  of  the  scaling.  From  FeriiaiKles  fl  al.  (1992). 

where  a  and  6  are  the  major  and  minor  axes  of  the  ellipse;  and  w  is  the  width  of  a 
wavenund>er  bin  in  the  2-d  s|>ertruni.  The  two-dimensional  spectrum  P{k,v)  for  the 
full'field  dopplergrams  is  shown  in  figure  (i.l. 
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6.2  Parametrization  and  Least-Squares  Fitting 


6.2.1  Fitting  Procedure 


The  frequencies,  widths,  and  heights  of  the  ridges  were  measured  as  functions  of 
wavenumber.  To  do  this,  constant  wavenumber  slices  of  the  power  spectra  were  indi¬ 
vidually  fitted  by  least-sc|uares  method.  Examples  of  slices  at  constant  wavenumber 
appear  in  figure  6.5.  A  single  isolated  mode  should  have  a  Lorentzian  profile. 

where  Pq  is  the  mode  power;  //q  is  the  mode  frequency;  and  F  is  the  fwhm.  The  high 
frequency  mock  modes  described  l>y  Kumar  et  al.  (1990)  should  have  a  sinusoidal 
profile. 


P(u)  « 


sin*(u;a:o/c) 


'  V";  t  1  •  2/  /  \ 

sm  (wa/c) 

where  Xq  is  the  source  dei)th;  is  the  acoustic  cutoff;  and  o  is  the  height  of  the  upper 
reflecting  surface.  Solar  processes  such  as  horizontal  flow  fields,  and  measurement 
effects  such  as  seeing  and  apodization  tend  to  smear  the  shape  to  produce  a  Gaussian 
profile. 


where  Pq  is  the  mode  power;  i/q  is  the  mode  frequency;  and  er  =  f'whm/2.36. 

The  ridges  lie  on  top  of  a  background  which  decreases  monotonically  with  wave- 
number.  In  the  photosphere,  the  granulation  is  caused  by  overshoot  of  convective  cells 
into  the  overlying  convert  ively  stable  layers.  The  power,  in  the  photosphere,  below 
the  /-modes  is  caused  by  this  convective  overshoot.  Extrapolation  of  calculations  of 
Nesis  &  Mattig  (1989)  to  the  low  chromosphere  suggests  a  similar  explanation  for  the 
background  power  there. 

Because  of  the  uncertaint}'  of  the  physical  description  for  the  background,  there 
is  no  proper  way  to  parametrize  it.  Nevertheless,  the  background  must  be  accounted 
for  to  avoid  contaminating  the  ridge  fits.  It  can  either  be  included  as  an  extra  term  in 
the  fits,  or  fit  independently  and  subtracted  from  the  power  spectrum  before  fitting 
the  ridges. 
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Figure  6.5:  Constant  wavenumber  slices  of  the  3-d  power  spectrum.  Data  are  plotted 
as  points  with  error  bars.  Background  and  individual  Gaussian  modes  are  plotted 
as  dashed  lines;  and  the  sum  of  background  and  modes  is  plotted  as  a  solid  line. 
Top  left:  k  =  1.16  Mm~';  Bottom  left:  k  =  5.25  Mm~*;  Right  (top  &  bottom): 
k  =  2.34  Mm-*. 
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The  latter  method  was  chosen  because  the  background  could  be  well  characterized 
using  the  range  of  frequencies  on  either  side  of  the  ridges.  The  shape  of  the  background 
P/s  suggested  a  power  law  paranietrization. 

PkM  =  (6.6) 

where  /Ifc  and  Bk  are  parameters  fitted  independently  for  each  wavenumber  slice. 
The  frequency  ranges  used  for  the  background  fits  include:  (0.12  <  v/mHz  < 
1.3)  and  (12.0  <  i//mHz  <  l.j.O). 

The  ridges  were  fitted  after  subtracting  the  background.  As  can  be  seen  from 
figure  6.5,  the  ridges  overlap,  and  so  must  be  fitted  simultaneously.  With  three 
parameters  per  Gaussian,  and  up  to  eight  visible  ridges,  the  fits  involved  up  to  24  pa¬ 
rameters  at  a  time.  This  large  parameter  space  combined  with  noisy  data  necessitated 
the  development  of  a  ver\  roluist  fitting  algorithm. 

The  fitting  program  is  described  in  the  next  section.  It  provides  a  nonlinear, 
weighted,  least-squares  fil  to  ungridded  data.  The  weights  used  in  the  least-squares 
fit  were  the  measurement  errors.  Hence,  the  fitting  procedure  minimized  x*-  The 
error  on  the  power  was  estimated  from  ecpiation  2..38.  This  in  turn  assumed  the 
velocity  and  intensit}-  errors  deri\ed  in  section  2.1. 


6.2.2  Levenberg-Marquardt  Method  for  Nonlinear  Regres¬ 
sion 


The  basic  algorithm  follows  the  Levenberg-Marquardt  procedure  (Press  et  al.  1989). 
The  method  attempts  to  minimize  the  in  parameter  space  (a),  given  N  measured 
values  y,  ±  <t,  at  points  .r,. 


(6.7) 


where  //(.r;rt)  is  the  fitted  function. 

The  minimum  of  the  surface  \*(n)  is  searched  for  in  an  iterative  manner.  Succes¬ 
sive  guesses  for  the  parameter  vector  a  are  made  based  on  the  gradient  and  curvature 
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matrix  for  the  suifac*'.  When  the  distance  to  the  minimum,  as  estimated  from 
the  curvature,  is  large,  the  gradient  is  weiglited  more  heavily  in  the  guess.  When  the 
minimum  is  near,  the  curvature  carries  more  weight.  The  algorithm  assumes  it  has 
found  the  minimum  when  it  can  no  longer  reduce  by  more  than  0.1.  The  covariance 
matrix  for  the  fit  is  simply  the  inverse  of  the  curvature  matrix  at  the  minimum. 

To  prevent  runawaj's.  the  parameters  are  tested  after  each  iteration  of  the  algo¬ 
rithm  in  which  is  reduced.  Parameters  that  are  not  within  set  limits  are  fixed  at 
appropriate  values,  and  no  longer  varied  in  the  fit.  The  remaining  parameters  are 
reset  to  their  initial  values  and  the  procedure  starts  again  from  the  beginning. 

For  Gaussian  ridges,  the  limits  for  the  height  parameter  are  e  (small,  positive) 
and  /i,„ax  (very  large).  X'alues  falling  outside  the  limits  are  fixed  at  the  limit.  The 
limits  for  the  width  parameter  are  one-half  and  twice  the  expected  width.  Again, 
values  falling  outside  the  limits  are  fi.xed  at  the  limit.  The  limits  for  the  frequency 
parameter  fall  halfway  between  the  expected  positions  for  ridges  (So).  For  the  bottom 
(top)  ridge  in  a  fit,  the  lower  (upper)  limit  is  defined  such  that  the  two  limits  are 
equidistant  from  the  expected  ridge  position.  Values  falling  outside  the  limits  are 
fixed  at  the  initial  guess  (Tq. 

The  reason  for  this  ticatment  of  parameter  limits  is  to  prevent  noisy  ridges  from 
disturbing  the  fits  of  adjacent  ridges.  Ridges  lohich  are  too  noisy  to  be  properly  fitted 
are  treated  as  background. 

For  the  first  spectrum  to  be  fitted,  initial  guesses  at  a  number  of  wavenumbers  were 
made  by  eye,  using  cross  hairs  on  a  workstation  display.  Guesses  at  other  wavenum¬ 
bers  were  interpolated  or  <'xtrapolated  from  the  nearest  four  fitted  wavenumber  slices. 
For  the  remaining  power  spectra,  initial  guesses  were  interpolated  from  the  already 
fitted  spectra. 

In  addition  to  the  background  and  ridge  fits,  the  least-squares  fitting  program  was 
used  to  interpolate  and  extrapolate  mode  frequencies  to  make  the  initial  guesses,  and 
to  determine  the  parameters  of  equation  6.1  in  section  6.1.1. 


k  (Mm"') 

Figure  6.6;  Povver  law  <^\pouent  for  background  fit.  Exponent  vs.  wavenumber. 

6.3  Results  of  Fits 

6.3.1  Background 

The  exponent  to  the  power  law  fit  is  plotted  in  figure  6.6.  The  background  velocity 
power,  integrated  over  frequency,  is  plotted  in  figure  6.7.  The  spectrum  falls  off  much 
more  slowly  with  wavenumber  than  the  spectrum  of  convection  in  the  photosphere  re¬ 
ported  by  Chou  et  ai  ( 1991 ).  The  root- mean-square  (RMS)  velocity  is  about  820  m/s. 
Nesis  k  Mattig  (1989)  have  measured  RMS  velocities,  due  to  ‘^mechanical  overshoot”, 
that  increase  with  height  in  the  upper  photosphere.  Their  model  predicts  the  height 
dependence,  for  Une-of-sight  velocities,  shown  in  figure  6.8.  A  smooth  extrapolation 
of  this  curve  to  the  h«ght  of  formation  of  the  Na  Dj  line  (500-700  km)  yields  an 
RMS  velocity  of  700-1100  m/s,  in  agreement  with  our  measurements. 


10®  10' 

k  (Mm"') 

Figure  6.7:  Background  velocity  power  integrated  over  frequency.  Background  was 
fitted  from  the  ‘‘quiet’’  siibfield  spectrum.  From  Fernandes  et  al.  (1992). 

6.3.2  Mode  Frequencies 

Frequencies  of  ridges  from  the  fuU-field  velocity  and  intensity  spectra,  and  the  quiet 
subfield  velocity  spectrum  are  tabulated  with  their  respective  statistical  errors  in 
appendix  A.  The  velocity  points  are  plotted  in  figure  6.9.  All  tabulated  points 
correspond  to  fits  in  which  the  offset,  width,  and  height  of  the  Gaussian  converged 
to  physically  reasonable  values,  and  the  integrated  power  under  the  Gaussian  was  at 
least  4  standard  deviations  above  the  background. 

All  six  fitted  p-mode  ridges  extend  up  to  at  least  7.0  mHz.  The  ps  and  ps  ridges 
extend  up  to  a  frequency  of  9.5  mHz.  There  is  no  deviation  from  a  smooth  extrap¬ 
olation  of  the  ridges  at  lower  frequencies.  Libbrecht,  Woodard  &  Kaufinoan  (1990) 
have  made  measurements  that  overlap  with  ours  at  frequencies  up  to  5.3  mHz.  As 
can  be  seen  from  figure  6.10,  the  fitted  frequencies  are  consistently  higher  than  those 
of  Libbrecht  et  al.  by  an  average  of  (35  ±  13)  pHz.  This  difference,  however,  is  only 
a  quarter  the  frequency  bin  width,  and  is  indicative  of  the  S3r8teinatic  error  in  the 
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Figure  6.8:  Calculated  Doppler  velocity  vs.  height  for  the  Fe  I  (5123.75  A)  line.  From 
Nesis  &:  Mattig  (1989). 

measurement. 

The  /-ridge  is  visible  at  wavenumbers  up  to  5.57  Mm"*.  There  is,  however,  a 
systematic  discrepancy  lietween  frequencies  fitted  from  the  velocity  and  intensity 
spectra  as  can  be  seen  from  figure  6.11.  Frequencies  measured  from  the  velocity 
spectra  are  consistently  higher  than  those  from  the  intensity  spectra  by  up  to  50  /iHz 
at  wavenumbers  below  2.  )  Mm"*.  This  difference  increases  to  ~  100  fiHz  at  4  Mm"*. 
The  discrepancy  may  in  part  be  due  to  a  difference  between  background  shapes  in 
the  velocity  and  intensity  spectra. 

Libbrecht  et  al.  (1990)  report  a  Gaussian  background,  centered  near  3  mHz,  in 
their  velocity  spectrum.  Inclusion  of  such  a  background  in  our  velocity  fit  would  bring 
the  /-ridge  frequencies  into  better  agreement  with  the  intensity  data  at  wavenumbers 
below  1.1  Mm"*,  but  would  push  them  farther  apart  at  higher  wavenumbers.  Fitting 
such  a  background  function,  moreover,  would  require  a  longer  time  series  with  corre- 
spcmdingly  better  signal- to- noise  ratio.  It  is  noteworthy  that  the  /-mode  frequencies 
measured  from  the  intensity  spectra  differ  with  those  rep<Mrted  by  Libbredit  et  a/,  by 
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Figure  6.9:  Measured  mode  frequencies  (from  velocity  spectra):  v  vs.  logtoib.  Fre¬ 
quencies  of  the /and  pi-pi;  ridges  are  plotted  with  1  o  error  bars.  For  k  <  3.82  Mm~*, 
frequencies  are  taken  from  the  full-field  spectrum.  At  higher  k,  values  (roai  the  quiet 
subfteld  are  plotted. 


CHAPTER  6.  MEASUREMENTS  FROM  POWER  SPECTRA 


87 


Figure  6.11:  Comparison  of  /-mode  frequencies  measured  from  velocity  and  intensity. 
Frequency  difference  vs. 

under  30  fiHz,  as  shown  in  figure  6.12. 

The  weighting  of  magnetic  regions  is  different  for  intensity  and  velocity  spectra. 
However,  as  shown  in  section  6.3.5,  the  discrepancy  is  not  caused  by  a  difference  in 
mode  frequencies  in  the  magnetic  and  nonmagnetic  regions  within  the  fidd  of  view. 
Regardless  of  which  set  of  /-mode  frequencies  we  use,  the  measured  frequencies  lie 
below  the  theoretical  frequencies  given  by  equation  1.9  at  wavenumbers  greater  than 
1  Mm~'.  The  difference  is  100-200  /i  Hz  at  4  Mm~’.  The  fractional  difference  between 
measured  and  theoretical  frequencies  is  plotted  in  figure  6.13. 

6.3.3  Mod«  Widths 

Ridge  widths  are  plotted  as  error  bars  to  the  frequencies  in  figure  6.14.  The  ridge 
widths  at  low  frequenqr  are  due  to  the  apodizing  function.  At  h%h  frequency,  the 
rMges  ate  broader  than  expected  for  mode  inodes.  The  measured  widths  ptobaUy  do 
not  represent  the  true  widths.  The  width  is  the  most  sensitive  parameter  to  mnse  in 
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Figure  6.12:  Comparison  of  /-mode  frequencies  from  full-field  intensity  spectrum  with 
Libbrecht  et  al.  (1990). 

the  fitted  power  spectrum. 

6.3.4  Mode  Power 

Mode  power  is  defined  as  the  integral  under  the  fitted  Gaussian. 

where  h  is  the  height  of  the  ridge;  and  w,  the  /whm.  The  error  on  the  power  is  given 

where  is  the  covariance  of  the  hdght  and  width.  The  mode  power  is  pk^ted 
as  a  function  of  frequency  in  figure  6.15.  The  p-mode  power  has  the  same  frequency 
dependence  for  all  p-ridges.  There  is  very  little  explicit  wavenuinl>er  dependence. 
Power  increases  as  to  a  maximum  at  3.4  mHz,  at  whi<&  p<nnt  it  turns  over  and 
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Figure  6.13:  Comparison  of  measured  /-mode  frequencies  with  theory.  Frequencies 
from  the  full  field  intensity  spectrum  are  compared  with  the  relation  w*  =  gk  {g  = 
274.1  m/a’). 
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Figure  6.14:  Frequency  and  ridge  width  v$.  k.  Top  plot  shows  /-modes  and  even 
order  p-modes.  Bottom  plot  shows  odd  ordered  p-taodes.  Error  bars  represent  the 
fwhtn  from  the  fit. 


CHAPTER  6.  MEASUREMENTS  FROM  POWER  SPECTRA  91 


t/  (mHz) 


Figure  6.15:  Velocity  power  vs.  mode  frequency.  From  Fernandes  et  al.  (1992). 

falls  off  again  as  The  /-modes  have  less  power  at  a  given  frequency  than  the 
p-modes.  As  surface  gravity  waves  with  very  little  compression  they  are  ph3^ically 
different  from  the  p-modes  and  are  not  expected  to  show  the  same  power  spectrum. 

An  MTF  corrwtion  was  made  to  the  power  plotted  in  figures  6.7  and  6.15  by  com¬ 
bining  the  calculated  MTF  of  the  telescope  and  CCD  with  the  estimated  atmospheric 
MTF.  See  sections  3.4  and  5.4.1. 

6.3.5  Comparison  of  Magnetic  and  Nonmagnetic  Regions 

Mode  frequencies  were  compared  for  the  quiet,  canopy,  and  sunspot  subareas  shown 
in  figure  6.1.  The  frequency  differences  are  plotted  in  figures  6.16  and  6.17.  The 
frequencies  of  ridges  in  each  subarea  agree  with  those  of  the  other  subareM  to  within 
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Figure  6.16:  Mode  frequency  differences  for  canopy  and  quiet  subfields.  Top  left, 
/-mode;  top  right,  pi;  bottom  left,  pa;  bottom  right,  pg. 


the  scatter  of  the  data  (±40  pHz). 

The  experimental  scatter  is  much  larger  than  the  error  bars  derived  from  the  fit. 
This  may  simply  indicate  an  underestimation  of  the  velocity  error,  or  it  may  be  a 
result  of  solar  noise.  Modes  of  finite  lifetime  have  an  intrinsic  frequency  uncertainty. 
Determination  of  the  frecjuency  of  the  cavity  resonance  i/©  requires  observations  for 
several  mode  lifetimes.  The  2.4  hr  observing  run  probably  constitutes  less  than  one 
mode  lifetime  over  much  of  the  wavenumbw  range.  Hence,  the  measured  frequencies 
will  be  distributed  about  i/©  in  a  manner  which  depends  on  the  mode  lifetime,  and 
shape  of  the  wavepackets. 
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Chapter  7 

Spatial  Distribution  of  Power 


The  spatial  distiiljiitioii  ol  power  was  measured  for  groups  of  modes,  as  well  eis  for  the 
sub-fundamental  region  of  tii<*  power  spectrum.  To  do  this,  the  set  of  dopplergrams 
V{x,t)  was  masked  using  ir(.r)  from  equation  2.33.  and  a  .3-d  FFT,  calculated  for  the 
full  field  of  view.  .A  sjiectial  mask  M{k.i/)  was  then  applied  to  the  Fourier  transform 
to  remove  all  but  a  small  range  of  wavenumbers  of  a  given  ridge.  An  inverse  Fourier 
transform  of  this  masked  3-d  spectrum  yielded  velocity  V\/(x,  <)  for  a  specific  set  of 
modes. 

Vxiir.l)  =  JF-' {.!/(?. //)jr3{lF(F)  1/(1,  <)}}  (7.1) 

The  mean  power  was  cah  ulated  at  each  pixel  to  make  a  velocity  power  map  Pm(x). 

Pm(^)  =  (7.2) 

=  j.  jy  ljF^'{3/(A^,,/)7^{lF(.f)V'(f,t)}}f  du  (7.3) 

where  <lenotes  the  spatial,  inverse  Fourier  transform. 

7.1  Mode  Masks 


Masks  were  constructed  fi»r  the  /  through  pa-modes  in  6  wavenumber  ranges.  In  addi¬ 
tion,  two  sub-fundamental  and  a  high  frequency  set  of  masks  were  made  for  the  same 
wavenumber  ranges.  Th<’  ma.sk  limits  are  shown  in  figure  7.1.  Three-dimensional 
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masks  were  made  from  these  2-d  limits  by  rotating  the  loci  in  the  azimuthal  direc¬ 
tion.  To  account  for  the  foreshortening,  the  rotation  followed  ellipses  described  by 
equation  6.1. 


7.2  Results 

Two  sample  power  maps  are  compared  with  a  magnetogram  and  an  H  a  image  in 
figure  7.2.  Only  the  central  disks  of  the  frames  are  visible  due  to  the  apodizing 
window.  The  subi-/  and  /-mode  maps  were  both  made  for  the  wavenumber  range 
1.95  <  A'/Mm”*  <  2..'M.  This  is  equivalent  to  a  wavelength  of  about  4",  or  3  Mm. 
The  mottled  appearance  of  the  power  maps  is  an  artifact  of  random  interference  of 
waves  in  a  finite  time  series.  The  subi-/  map  shows  little  variation  other  than  an 
overall  trend  of  higher  power  in  the  upper  right  hand  corner.  The  /-mode  map  shows 
a  depression  in  power  in  the  active  region. 

A  selection  of  power  maps  is  shown  in  figure  7.3.  An  estimate  of  the  RMS  velocity 
in  each  map  can  be  obtained  from  figure  7.4.  The  power  in  high  wavenumber  modes 
is  depressed  in  the  active  region.  .-\l  low  wavenumbers,  the  active  region  has  little 
effect  on  the  power.  1  he  sub|-/and  subi*/  niaps  show  no  correlation  of  power  with 
magnetic  field  at  any  wavenumber.  The  high  frequency  maps,  though  not  shown  in 
figure  7.3,  looked  very  similar  to  p-mode  maps  of  the  same  wavenumber  range  at 
lower  frequency.  There  i'^  little  difference  between  the  /,  pi,  p2,  and  pa  maps  for  a 
given  wavenumber  range. 

A  horizontal  slice  through  the  /-mode  map  from  figure  7.2  was  made  above  the 
sunspot.  The  slice  was  low  pass  filtered  using  a  Gaussian  impulse  response  function 
(16  bin  fwlim).  See  figure  7.').  This  low  pass  filtered  slice  is  compared  with  the 
magnetogram  signal  in  figure  7.0.  The  power  depression  is  significantly  broader  than 
the  spatial  extent  of  the  photospheric  magnetic  field.  It  is  narrower  than  the  canopy 
field  at  ~  1200  km.  as  can  be  .seen  from  the  H  a  image  in  figure  7.2.  The  filaments 
visible  in  H  o  are  streams  of  gas  traveling  along  magnetic  field  lines.  The  filaments, 
thus,  trace  out  the  field  lines  which  form  a  canopy  above  the  photosphere  and  lower 
diromosphere. 
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Figure  7.1:  Masks  for  power  distribution  maps.  Measiired  mode  frequencies  are 
shown  as  dots.  Mask  limits  are  shown  as  solid  lines.  The  high  firequency  (HF)  masks 
covered  the  range  6.0  <  <//mHz  <  10.0;  the  p-mode  mask  limits  were  set  halfway 
between  nunles;  the  /-mode  limits  were  set  equidistant  from  the  /-mode;  the  sub]-/ 
band  extends  from  the  /-mode  down  to  v  ~0;  the  8ubi-/band  is  the  ssmie  width  as 
the  /-mode  band.  It  is  a  subset  of  the  sub]-/ band. 
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Figure  i.2:  Com|)arisoii  of  |)o\vci  maps  witli  magnetic  fields.  Top  left;  log  power 
for  sul>i-/ (1.9.)  <  k/Wwi  '  <  2.;M).  Top  right:  absolute  value  of  magnetogram  (log 
scale).  Bottom  left;  Ho-  100  m.A  image.  Bottom  right;  log  power  for /-mode 
(1.9-5  <  A’/Mm  ‘  <  ‘i.-M).  1  he  small  tick  marks  are  5''. 
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Figure  7.3;  Power  rlistriliutiuu  for  select erl  masks.  Log  power  is  displayed  for  sub2-/, 
fi  Pi*  P2i  Pj  niaps  willi  wavemiinher  ranges; 

(1.95,2.34).  (1.17.1.56).  [(1.78.1.17).  (0.51.0.78)  Mm"'. 
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Figure  7.4:  RMS  velocity  vs.  k.  Curves  fitted  from  quiet  subfield  velocity  spectrum. 
Curve  labeled  sub-/ is  taken  from  the  background  fit. 

It  is  surprising  that  the  sub-fundamental  power  is  not  correlated  with  magnetic 
field.  The  RMS  convective  power  in  the  photosphere  is  a  factor  of  2  lower  in  plage 
(network  field)  than  in  nearby  quiet  regions  (Title  et  al.  1992).  Moreover,  the  subj-/ 
band  lies  just  below  the  /-mode  band,  but  shows  completely  different  behavior.  At 
the  same  wavenumber,  and  almost  the  same  frequency,  oscillation  amplitudes  are 
reduced  in  active  regions,  whereas  non-oscillatory  motkm  is  not. 

It  is  impossible  to  distinguish  between  instrumental  backgrounds  and  solar  back¬ 
grounds  from  a  single  data  set.  It  is  difficult,  however,  to  come  up  with  an  obser¬ 
vational  effect  that  would  produce  the  measured  badcground.  Filter  transmisskm 
variation  should  affect  only  the  very  low  wavenumber  part  of  the  spectrum.  Seeing 
variations  should  be  independent  of  frequency;  whereas,  the  observed  badcground 
decreases  with  frequency. 
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Chapter  8 

Conclusions  and  Discussion 


8.1  Surfaces  and  Cavities 

The  existence  of  a  chroniosplioric  cavitj’  is  expected  to  produce  a  resonance  consist¬ 
ing  of  a  horizontal  baud  of  jjower  and  -avoided  crossings”  at  about  5.5  mHz.  (See 
figure  1.4.)  This  cavity  is  reiiuired  by  the  model  of  Balmforth  &  Gough  (1990)  to 
explain  the  existence  of  higli  frc<iuency  inodes.  The  p-mode  and  /-mode  ridges  shown 
in  figures  6.4  and  6.9  are  <  learly  visible  at  5.5  mHz.  but  there  is  not  a  horizontal  band 
of  power  nor  are  there  "avoitled  crossings".  The  clear  absence  of  this  resonance  in  our 
data  shows  that  there  is  no  chromospheric  cavity.  Either  the  reflection  coefficient  at 
the  base  of  the  corona  is  too  low  or  the  surface  is  too  irregular  to  form  a  good  cavity. 

This  result  supports  the  alternate  hypothesis  of  Kumar  et  al.  (1990),  that  the 
high  frequency  ridges  are  a  purely  geometrical  artifact.  The  result  is  in  agreement 
with  the  conclusion  of  Woodard  Libbrecht  (1991);  and  it  contradicts  the  cl2dm  of 
Dame  et  ai  (1984)  that  there  is  a  horizontal  ridge  of  power  at  5.5  mHz.  Our  power 
spectrum  has  a  significant  ly  better  signal  to  noise  ratio  than  that  of  Dame  et  ai 

The  /-modes  are  sensitive  to  conditions  in  the  chromosphere  and  transition  re¬ 
gion.  The  fact  that  /-mode  frecpiencies  are  lowered  at  high  wavenumber  may  indicate 
interaction  with  chromospheric  magnetic  fields,  energy  leakage  into  the  corona,  or 
horizontal  inhomogeneiti«-s  in  the  chromosphere-cortma  boundary.  The  last  of  these 
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three  possibilities  is  interesting  from  the  standpoint  of  acoustic  cavities.  Perhaps  ir¬ 
regularities  in  the  interfai  e  lictween  the  chromosphere  and  corona  are  responsible  for 
both  lowering  the  /-mode  fre((acncies,  and  for  reducing  coherent  reflection  of  acoustic 
waves. 


8.2  Waves  in  Magnetic  Fields 

There  is  no  freciuency  shift  for  /  and  p-modes  in  magnetic  network  fields  of  100  — 
—1000  gauss  to  within  tin*  scatter  of  the  data  (~  50  pHz).  This  limit  is  not  stringent 
enough  to  determine  whether  tlie  freciuency  shifts  measured  by  Woodard  et  al.  (1991) 
result  from  small  shifts  over  the  whole  disk,  or  larger  local  shifts,  in  active  regions, 
which  contaminate  the  full-disk  measurement.  The  modes  measured  by  Woodard  et 
al.  have  widths  of  order  1  //Hz.  Hence,  to  contaminate  such  a  measurement,  mode 
frequencies  in  active  regions  would  have  to  be  distributed  within  a  few  /iHz  of  the 
quiet  sun  values.  .A  measurement  of  such  precision  would  require  observations  of  at 
least  20  hr  duration,  assuming  the  waves  are  coherent  for  the  observing  period.  .At 
La  Palma,  20  hours  of  observations  can  be  made  in  two  days.  Hence,  the  waves  would 
add  coherently  if  the  mo<le  lifetime  is  >  34  hr.  If  the  mode  lifetime  is  less  than  this, 
the  signal  to  noise  ratio  increases  only  as  the  scpiare  root  of  time,  so  a  much  longer 
observing  time  would  Ire  necessary. 

Power  maps  for  oscillations  in  the  low  chromosphere  show  that,  at  high  wavenum¬ 
ber,  wave  amplitudes  ar<'  suppressed  15  Mm  outside  of  the  photospheric  magnetic 
field.  (See  figure  7.6.)  In  the  photosphere,  the  power  is  suppressed  within  one  (hori¬ 
zontal)  wavelength  from  the  field  (Tarbell  et  al.  1988).  The  vertical  distance  between 
the  photospheric  and  chromospheric  observations  is  under  0.5  Mm,  and  the  (true) 
wavelength  of  the  oscillat  ions  is  2  Mm.  Diffraction  cwnot  account  for  the  spread¬ 
ing  of  the  region  of  the  power  depression  because  the  waves  have  traveled  less  than 
one  wavelength  in  the  interval  over  which  the  spreading  occurs.  The  power  suppres¬ 
sion  must,  therefore,  be  produced  by  interaction  of  the  oscillations  with  the  field  in 
the  layer  in  where  the  warts  an  obserrtd. 

The  area  containing  magnetic  flti.x  was  not  measured  directly  at  the  height  of  the 
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observations.  Its  extent  ran  lie  inlerred.  however,  from  the  magnetogram  made  in 
the  mid-photosphere,  an»l  the  H  a  image  in  the  mid-chromosphere.  The  field  lines 
spread  out  rapidly  above  the  photosphere.  Hence,  the  extent  of  the  field  at  the 
observed  height  should  be  intermetliate  between  the  limits  of  the  mid-photospheric 
and  mid-chromospheric  fields.  This  is  consistent  with  the  above  conclusion. 

The  /-modes  propagate  without  compression,  so  they  produce  no  pressure  fluctu¬ 
ations  in  the  region  where  oscillations  are  suppressed.  The  mechanism  suppressing 
the  wave  motion  must,  therefore,  act  through  the  velocity  fluctuations. 


8.3  Directions  for  Future  Work 

The  frequency  measiireni<  nts  presented  here  show  the  existence  of  high  frequency  and 
high  wavenumber  oscillations.  They  show  general  characteristics  of  the  dispersion 
relations  such  as  the  smooth  cotttinuation  above  the  acoustic  cutoff.  They  are  not, 
however,  high  precision  iiieasureinenls  in  comparison  with  existing  low  wavenumber 
observations.  Precision  of  the  measurements  can  be  improved  with  longer  observing 
times.  This  would  improve  the  signal  to  noise  ratio,  and  enable  better  fitting  of  the 
background  and  ridges.  With  better  frequency  resolution,  it  would  become  necessary 
to  correct  the  fitted  frequencies  for  their  weighting  with  mode  power  (Korzennik 
1991). 

Important  characterisi  ics  of  the  solar  oscillations  can  be  inferred  from  their  phases 
(Deubner  &:  Fleck  1989.  1990;  Fleck  &  Deubner  1989).  Trapped  and  nontrapped 
oscillations  can  be  distinguished  l)y  the  phase  difference  between  velocity  and  intensity 
measurements.  Propagating  and  evanescent  waves  can  be  distinguished  by  the  phase 
lag  between  velocities  at  (lifferent  heights  in  the  atmosphere.  An  attempt  was  made 
to  measure  velocity-intensity  phase  differences.  The  attempt  failed,  possibly  because 
of  the  large  height  range  included  in  the  measurement  (see  section  1.4),  uid  possibly 
because  of  the  large  number  of  modes  included  in  a  single  frequency  bin  (Braun  et 
al.  1992a). 

The  interaction  of  magnetic  fields  with  oscillations  can  be  better  quantified  by 
measuring  both  wave  amplitudes  and  magnetic  fields  at  a  range  of  heights  in  the 
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atmosphere.  Improvements  that  liave  been  made  to  the  tunable  filter  system  since 
the  1990  observing  season  will  allow  Doppler  measurements  i-estricted  to  either  the 
core  or  wings  of  the  Na  Di  line.  In  addition,  magnetograms  can  be  made  in  the  core 
or  wings  of  the  line  by  finding  the  difference  between  dopplergrams  made  with  left 
and  right  circularly  polarized  light. 

Braun  et  al.  (1992b)  report  halos  of  enhanced  high  frequency  acoustic  power 
around  active  regions.  \  search  was  made  for  such  halos  in  power  maps  made  with 
the  HF  masks  shown  in  figure  7.1.  No  enhancement  of  power  was  found.  This  does 
not  contradict  the  result  of  Braun  e!  al..  as  they  measured  intensity  fluctuations  in 
the  photosphere,  in  a  lower  wavenumber  regime.  It  would  be  interesting  to  repeat 
the  analysis  with  continuum  brightness  measurements  for  a  field  of  view  of  at  least 
6  arc  minutes. 

The  cause  of  the  sub-fundamental  background  is,  currently,  unknown.  Experi¬ 
ments  should  be  done  to  disc  riminate  between  observational  and  solar  backgrounds. 
For  example,  a  set  of  filtergrams  taken  at  fixed  wavelength  with  a  high  cadence 
could  be  used  to  make  anificial  dopplergrams,  and  check  the  effect  of  seeing  on  back¬ 
ground.  Another  set  of  filiergrams  taken  at  the  nominal,  observing  wavelengths  could 
be  made,  using  integrated  sunlight,  to  check  the  effect  of  filter-transmission  variation 
on  background. 

An  attempt  was  made  to  find  the  oscillation  sources  described  by  Brown  (1991).  A 
temporally  high-pass  filt<'red  movie  of  velocities  was  scanned  visually  for  evidence  of 
spatially  and  temporally  localized  sources  of  power.  No  good  candidates  were  found. 
A  better  method  might  b<*  to  construct  the  two-dimensional  acoustic  Pointing  vector 
described  by  Bravtn  et  al.  ( 1992b).  Large  positive  divergences  of  this  vector  field  for 
high-pass  filtered  data  could  point  out  oscillation  sources. 
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Table  A.l:  Mode  frequency  u  ±(Ju  (^Hz)  as  a  function  of  wavenumber  k  (Mm“^)  — 
measured  from  dopplergrams.  For  k  <  3.82  Mm~\  frequencies  are  taken  from  the 
full-field  spectrum.  .4t  higher  k,  values  from  the  quiet-subfield  are  tabulated. 
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Table  A.2:  /-mode  frequency  ±  (m^z)  as  a  function  of  wavenumber  k  (Mm“‘) 

—  measured  from  intensities. 
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